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ABSTRACT 


A  capillary  viscometer  was  developed  to  investigate 
the  rheological  behaviour  of  sand-water  suspensions. 

The  suspension  was  found  to  behave  as  a  non- 
Newtonian  dilatant  fluid  over  the  range  of  shear  rates  tested. 
Bagnold 11  s  concept  of  grain  shear  stress  T  and  linear  con- 

centra t ion  A  are  used  in  dimensional  analysis  of  the  flow 

2-  2 
i  Tf  d  r  sv  \  d 

situation.  The  dimensionless  numbers  anc*  v_Dj~Yj—  are 

found  to  bear  unique  functional  relationship  for  all  values 

of7^<3.5.  By  comparison  with  the  experimental  observation 

2. 

of  Bagnold,  the  existence  of  a  critical  ^ is  postulated  such 

-2. 

that  the  G)  marks  the  transition  from  Newtonian  to  non- 

Tcv 


Newtonian  flow  regime. 

A  critical  concentration  at  which  there  is  an 
abrupt  change  in  the  characteristics  of  the  suspension  is 
noticed  at  concentrations  between  0.30  and  0.35. 
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INTRODUCTION 


1. ,  X  NON- NEWTONIAN  FI  OW 

In  recent  years  there  has  been  a  marked  increase 
in  interest  in  the  flow  and  related  behaviour  of  complex 
mixtures,  especially  of  solid- liquid  suspensions.  Hydraulic 
transport  of  solids  is  used  extensively  in  road  and  dam 
construction,  land  and  harbour  reclamation,  in  mining  and 
chemical  industries.  In  connection  with  the  research  in¬ 
vestigations  or  design  work  on  the  problems  of  hydraulic 
transport  a  knowledge  of  the  consistency  or  the  "viscosity" 
of  the  suspension  is  usually  called  for.  However,  most  of 
the  suspensions  happen  to  be  non- Newton i an  in  character. 

The  study  of  the  rheological  behaviour  of  non-Newtonian 
fluids  in  general  and  the  suspensions  in  particular  is  an 
extremely  complex  and  incompletely  understood  field.  Advances 
are  being  made  in  the  various  disciplines. 

The  object  of  the-  present  investigation  is  to 
study  certain  aspects  of  the  rheological  behaviour  of  sand- 
water  suspensions.  To  understand  the  present  status  of  the 
know ledge  in  this  topic,  it  is  necessary  to  know  the 
characteristic  behaviour  of  non- Newtonian  fluids.  A  brief 
review  of  non-Newtonian  fluid  flow  is  given  in  this  chapter. 

1 , 2  NON  -  NEWT, OK  I  AN  FUJI  D,S  : 

The  behaviour  of  fluids  when  subjected  to  shear 
stress  forms  the  basis  for  classifying  the  fluids  into 


Newtonian  and  non- Newtonian  categories. 


:i 


2 


Fluids  for  which  the  shear  rate  is  directly 
proportional  to  the  imposed  shear  stress  only  are  known 

as  Newtonian  fluids.  The  rheological  equation  is 

^ 


or 


T  -  "1  Sf  - CD 

Newtonian  behaviour  has  been  found  to  be  common  in  all 
gases  and  liquids  or  solutions  of  low  molecular  weight 


materials. 

All  fluids  which  do  not  obey  the  law  of  direct 
proportionality  of  T  and  ^  in  equation  1  are  non-Newtonian 
fluids.  They  are  further  classified  as: 

(a)  Time  independent  or  purely  viscous  fluids: 
Fluids  for  which  shear  rate  depends  only  upon  shear  stress. 


i 

(b)  Time  dependent  fluids: 


(2) 


Fluids  for  which  shear  rate  depends  upon  imposed  shear 
stress  and  the  duration  of  the  stress. 


i  -  (t  ,  t ) 


(3) 


(c)  Visco-Elastic  fluids: 


Materials  for  which  shear  rate  depends  upon  imposed 
stress  and  extent  of  deformation. 


i  -  C^E) 


(4) 


(d)  Complex  Rheological  bodies: 


Systems  exhibiting  the  complexities  of  the  categories  a»b. 


and  c  materials  in  combination. 


3 


(e)  Fluids  for  which  shear  rate  depends  upon 
shear  stress,  and  intensity  of  an  imposed  magnetic  field. 

(f)  Fluids  which  must  be  considered  as  being 
made  up  of  discrete  particles  rather  than  as  a  continuous 
me  d i a . 

Non-Newtonian  behaviour  is  commonly  observed  in 
(i )  solution  or  melts  of  high  molecular  weight  polymeric 
materials,  (il)  suspension  of  solids  in  liquids  at 
sufficiently  high  concentration  to  affect  or  associate 
with  the  liquid  phase. 

Even  though  there  are  six  categories  of  non- 
Newtonian  fluids,  virtually  all  engineering  work  to  date 
has  dealt  with  the  first  category.  Also  since  the  problem 
of  interest  in  the  present  investigation  is  the  one 
concerning  this  category  the  fluid  behavior  of  only  the 
Time- Independent  non-Newtonians  is  reviewed  here. 

1.2,1  TIME  INDEPENDENT  NON- NEWTONIAN  FLUIDS 

The  rheological  equation  'jf  =  holds  for  this 

category  of  fluids.  They  are  further  sub-dividad  into  3 
distinct  types  depending  on  the  plotted  form  of  the 
equation.  Figures  1  and  2  show  the  classification  by  the 
so  called  "flow  curves". 

(a)  Bingham  Plastic: 

This  is  the  simplest  of  non- Newton i an  fluids  and  has  the 
rheological  equation 

T-Ty  for  T>Ta  . (5) 

where  ^ ^  =  plastic  viscosity  or  coefficient  of  rigidity 


Bingham  Plastic 
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Fig.  1  Fluid  Flow  Curve  (Arithmetic  Plot) 


Shear  Rate  y 

Fig.  2  Fluid  Flow  Curve  {  Logarithmic  Plot  ) 
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-  initial  shear 
movement  i.e.  X  "  °  for 
The  flow  curve 


stress  necessary  to  initiate 

T  <  Tb 

is  shown  to  ordinary  and 


logarithmic  co-ordinates  in  Figures  1  and  2.  To 
Logarithmic  Co-ordinates  the  curve  becomes  asympotltic 
to  TJ  at  low  shear  rates  and  approaches  a  slope  of  unity 
at  high  shear  rates* 

The  explanation  for  Bingham  plastic  behaviour 
assumes  that  the  fluid  at  rest  contains  a  three  dimensional 
structure  of  rigidity  sufficiently  great  to  resist  finite 
stress  When  the  stress  is  exceeded  the  structure  breaks 

down  completely  and  the  fluid  displays  the  linear  relation- 

* 


ship  between  shear  stress  excess  and  the  shear  rate  (1)* 

This  concept  of  an  idealised  Bingham  plastic  is  very 
convenient  in  practice*  Even  though  there  is  no  convincing 
evidence  to  support  the  suggestion  that  any  real  fluids  obey 
the  Bingham  plastic  law. (2)  the  behaviour  of  many  real 
fluids  can  be  approximated  with  this  law*  Drilling  muds, 
chalk  suspension,  oil  paints,  toothpastes  and  sewage  sludges 
have  been  stated  to  follow  the  Bingham  law  approximately  (1*3'< 
(b)  Pseudo-Plastic i 

For  these  fluids  the  'Apparent  viscosity"  -  the  ratio  of 
shear  stress  to  shear  rate  -  falls  progressively  with  shear 
rate  and  reaches  a  limiting  value  at  very  high  rates  of  shear. 
Hence,  the  slope  of  the  logarithmic  slope  curve  is  less  than 
unity. 


Numbers  in  brackets  (  )  refer  to  the  list  of  references* 


O  -  v  ,  9  .  i 

>J  : 

-  - 

•  <  i  i  </  •  g  •  •  ' 

•  •  •  '  - 

.  .  )  n  ■  .  ■  '  •  ■  . 

'  ■ 1  i  •  <  .  •  ■  *< 

■ 

1  «  r:  o:. 

■■  "■  •  *  <  *'  .  ••  «r 
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(c)  Dilatant  fluids: 

The  apparent  viscosity  for  these  materials  increases  with 
increasing  rates  of  shear.  The  slope  of  the  logarithmic 
flow  curve  is  greater  than  unity. 

To  portray  the  behaviour  of  time  independent 
non-Newtonian  fluids  in  general  and  p seudo- p la s t i c s  in 
particular,  a  large  number  of  empirical  equations  have 
been  proposed.  A  very  popular  and  versatile  equation  which 
predicts  the  behaviour  approximately  over  a  range  of  shear 
rates  is  the  '’Power  Law" 

(6) 

where  K  and  n  are  constants  for  a  particular  fluid:  K  is 
generally  considered  as  a  measure  of  the  consistency  of  the 
fluid,  n  is  a  measure  of  the  degree  of  non-Newtonian  behaviour. 
This  equation  with  only  two  constants  is  used  for  pseudo¬ 
plastics  (n<l)  and  also  for  dilatant  fluids  (n>l). 

This  powerlaw  has  three  basic  defects,  as  pointed 
out  by  Reiner  (4).  Since  the  dimensions  of  T  and  ^  are 
fixed,  the  term  K  has  dimensions  which  depend  upon  the 
value  of  n.  It  is  assumed  that  different  materials  will 
have  different  value  of  n,  and  hence  K  is  different  both 
qualitatively  and  quantitatively  for  different  materials. 
Therefore  K  cannot  be  a  measure  of  a  single  property  such 
as  consistency. 

T_  . 

The  apparent  viscosity  Ick.~  ~  ^  ^  *  Hence 

>.  0  i.e.  when  fluid  is  at  rest, and  n<l. 


T  =  v<  C  V) 


^0  - ►  e>o  as  i 

'-a. 


■  J 


f 


This  is  contrary  to  experience  as  any  fluid  must  have  a 
definite  apparent  viscosity  at  rest;  otherwise  it  cannot 
begin  to  f low, 

Also  when  n  <  1  ,  ^  . — *-o  as  1$  — *-oo  wh ich  is  object¬ 
ionable  as  there  is  no  liquid  whose  viscosity  vanishes. 
These  objections  which  have  been  discussed  in 
detail  by  Reiner  (4)  are  assumed  not  to  be  serious  in  most 
engineering  applications  (3)  and  the  power  law  continues 
to  be  used  in  the  same  form  as  equation  6,  The  other 
empirical  equations  are  more  difficult  to  use  than 
fowerlaw  and  usually  do  not  offer  any  compensating 
advantage  s  ( 3 ) , 

FLOW  Op  NON-NEWTONIAN  FLUIDS  THROUGH  PIPES, 

The  case  of  steady,  laminar  flow  of  a  time- 
independent  non-Newtonian  fluid  through  a  circular  pipe 
is  discussed, 

Consider  a  fluid  for  which  is-'f(^)  .  For  flow 
through  circular  pipes  ^  «  By  considering  the 

equilibrium  of  a  cylindrical  element  as  shown  in  Figure  3, 


?,TT  V  L  T  =  TV  ^ 

_  Jf_ 

^  ~  L  2. 


At  the  vail 

T 

r 

L  9- 

T 

-  T  — 

'u*  R 

is 

■&U. 
”  ”  *5Y 

=  f  m 

=  fCT„ 

k) - 

Integrating, 

u(r) 

-  u(R )  - 

Tk  v)  A 

« 

3*  J  nr!  , 

■ 


8 
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If  there  is  no  slip  at  the  wall  u(R)  -  0 


v-00  -  |  (  C  -^r  Y  )  Ai 


Mean  velocity  V  •  11 R. 


(R 

\  ZTT  Y  U.(_Y)  cly 


U.  CY  )  Cl 


Uy1) 


I 

R> 


^  V2"  d  'u.OO 


But 

and 


u(R)  -  0 

d  =  f  CT-  -£-)  A, 


from  equn*  (7) 


t  [  _  i-  r  ,  • T  V  ,  _ 

w  -=r-  )  <=W 


V  ^ 


)  ^ 


(8) 


Again  Y  =  R 


T_ 

T, 


so  that 


A  y  =  —  clx 


T 

v*  i . 


Hence  inserting  in  equn*  (8) 


Lu> 


V  =  R 


t; 


T*  {(T)  C^T 


<8  v 


or 


*  \  T^fCO  cIt 


D  T 


(9) 


CO 


This  is  the  general  expression  for  the  relation¬ 
ship  between  the  "flow  function"  j§X  and  the  boundary  shear 
for  the  case  of  steady  laminar  flow  through  circular  pipes* 
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Knowing  the  form  of  function  f(^)  the  above 
equation  can  be  integrated  to  solve  for  various  pipe 
flow  parameters  in  terms  of  others* 


1.3*1  FOR  A  NEWTONIAN  FLUID 


i  •  e  • 


From  equation  (9) 

“5V 

"D" 


or 


But 


_D 

L  -4 


■52,  Y|  v  L 

5* 


(10) 


-  -  Hagen  - 

Poiseulle  equn.(lOa) 


Also  defining  the  friction  factor  f  as 


from  (10a) 


££  -  c  l  vz 


f 


64 
^  V  D 

n 


(id 


Where 


Re  = 


Reynolds  number 
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The  relation  Tw  oYlllj)  suggests  that  a  plot  of 
Tw  vs  <  )  can  be  used  for  testing  the  cases  other  than 

laminar  and  Newtonian  since  deviation  from  straight  line 
in  such  a  plot  indicates  a  breakdown  of  Newtonian  laminar 
conditions. 


la.  3, 2. 


FOR  A  BINGHAM  PLASTIC  FLUID, 

T  -  T-y  -  y  -  -  -  for  T  -  -  equn(5 ) 


'S  -  O 


-  -  -  for 


T  <  Tv 


if  ^  -  Tf  (,T‘Ta) - forT>t, 


n 

f  en  ■-  ° 

Equation  (9)  becomes 


f  or  T  <  T 


BV 

r> 


T 


Tj  -  v°° 

^  T,L  -f  (T)  +  (  T"  {C'O  <At 


w  o 


*03 


Tv 


BV 

X> 


T 


10 


T  CT-T.J  ^ 

X 


D 


4- 

-  tj*  4 

*-  UJ 

%  T 

< 

Y]  T° 

Ip  Vu> 

4- 

\ 

3 

\<z 

1 

-  ^  ( 

^ 

Vi!  \ 

- 

3  V 

Tuj  ) 

'  p 

(12) 


This  equation  (12),  known  as  the  Buckingham  equation, 
enables  the  relationship  between  pressure  drop  and  flow  rate 
to  be  determined  provided  the  fluid  properties  Ty  and  ^are 
known.  Note  that  it  reduces  to  equation  (10)  when  Ty**0,  and 
generalise  equation  (10)  to  allow  for  viscosity  Yj  being 
replaced  by  two  parameters  Ty  and  Y|  « 
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To  solve  for  the  pressure  loss,  the  Buckingham 
equation  has  been  expressed  in  non-dimensional  form  in  a 
number  of  ways,  Govier  and  Winning  (5)  divide  the  equation 

-  el— 

=&  ?  v 


(12)  by 


n, 


to  introduce  the  friction  factor  in  a  general¬ 


ised  equation  (11):  - 

l  {  Y 

Re'on  <c>4  <3 


64-  Y 
•3  f3  Rg 


where  f 

R 


T> 

L 

VDf 


t  V  /  2. 


"Ip 

Tt,  d 

vyIp 


-  4  (13) 

=  friction  factor 

-  a  modified  Reynolds 
number 

*  a  "yield  number" 


Perkins  and  Glick  (1.5)  have  written  equation  (13) 


a  s  : 


_i 

R 


erm, 


•f 

<84 


H 


G4  H 


6  R 


ev>^ 


■3)  f3  R 


B 


- (14) 


where 


H  “  - — -  »  the  Hcdstrom  number 


It  will  be  observed  that  11=  Y  K  .  Charts  are  available 

cm 

(1,3)  to  give  the  variation  of  f  with  H  (or  Y)  and  Rem  for 
the  laminar  region  -  as  given  by  equation  (12)  and  for  the 
turbulent  region  also  -  as  determined  by  experimental 
observa  t ions# 

1 , 3, 3  FLUIDS  WHICH  OBEY  POWER  LAW , 


If 


T  -  K 

T/ 


(6) 


x  -  (  7k)'"  =  Ccx) 

From  equation  (9) 


i  *  ®« 


bv 

4_ 

\ 

D 

T  3 

v 

4 

D 

KVn 

^  CO 

4m 

( 

T 


T 


2.-,  — 
n 


uj 


■+  '/n 


Tuj  \'/n 


3n  -+ » 


K 


) 


-  -  -  -  (15) 


13 


or 


T 


w 


(  A  , 

r 

- (16) 


or 


A|p  .  4^  (j 


r\ 


0 


^  . 

v  L 


D 


n-*-» 


( 16~a  ) 


It  is  interesting  to  observe  from  this  equation 


v"L 


(17)  that 


A\3  ^ 


i .  e 


D 

r> 


v\  -*■ 1 


2>n-H 


A  cc  L  j 

i.e,  For  a  given  rate  of  flow  in  a  known  length  of  pipe  the 
relation  between  the  diameter  of  the  pipe  and  the  correspond¬ 
ing  change  in  the  pressure  depend  upon  r>  of  the  fluid*  Thus 
for  pseudo-plastics  (n<D  the  pressure  drop  is  less  sensitive 
to  diameter  than  for  dilatant  fluids  (n>l).  This  fact  is 
important  in  the  design  of  capillary  viscometers  and 

non- 

metering  devices  f or^ Newtoni an  fluids. 

The  equation  (6)  together  with  equation  (15)  can 

~c>  \a,  . 

be  used  as  follows  to  get  an  expression  for  (  —  :r~  )  i*e» 


the  shear  rate  at  any  radial  distance  r 

/-  T  \'/n 


but 


i  -  -  -  C  V) 


fro©  eqim  (6) 


fe-c’sf'O) 


Vo 


Using  equation  (15) 


,  *2>vx  \  sn  +  l  / 

(r  ^  J  =  *4^“^  “5"K  *) 


Vo 


(17) 


1,3.4  ANY  TIME  INDEPENDENT  NON- NEWTONIAN  FLUID 


V  -  %  j  V  f(v  ^ 


(18) 


is  a  function  of  »  say  (Too) 

=  i,  (  f  -f(T)  dT 

*-00  J 

_  \  4[5f+CT«)] 


4  (T-) 


4  Tb 


cl  TuJ 


. 


■ 
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£«e«  f  ~  ^  ct^CTo) 

=  I  tc^)  +  5  44^ 

\  d  luj 

But  ■f  (.To)  -  iM  =  C-4L 

r_  _  2_  (  +  1 h  <H  %0 

^  _&Y  4  4  ^  s  ’  4  oL  tw  ,  , 

- (19) 

This  equation  which  gives  the  shear  rate  at  the  wall  for 
any  time  independent  non-Newtonian  is  known  as  the  Robin- 
owitsch  -  Mooney  equation.  The  solution  as  above  was  first 
given  by  Mooney  (6S1). 


Metzner  and  Reed  (7)  have  made  the  following 
useful  alteration  of  the  above  equation  (19)  using  the 


i dent i ty 

Tw 

A(^) 

_ 

1  4  0*  ¥) 

4 

d  Xui 

4 

4  Lo<g  Tea 

4" )  = 

v  ^  L 

3  /  8V 

'  4  V  X) 

)  +  4  ( 

QV', 

35  A  Lo^Tm 

Let 

d 

8V  \ 

) 

i 

A  ^ 

To 

r\' 

------  (20) 

that  i,« 

is  the  slope  at  any 

point  of 

the  curve  logT  vs 

Then  ^ 

"c)\A_  \ 

(  3n4l 

)  (t?) 

bY  ) 

7  00 

\  .4  n' 

------  (21) 

Also  from 

equation 

(20),  log 

Tw  -  n1 

log  Oj^)  +  Log  K^- 

T  ,  \  /  BV  > 

^  L  X)  J 


•  • 


-  -  (22) 


■ 


. 
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Equations  (21)  and  (22)  should  hold  good  for 

any  time  independent  fluid  regardless  of  the  shear  rate  - 

shear  stress  relationships  provided  there  is  no  slip  at 

the  wail.  Thuss  even  though  the  shear  stress  and  shear 

rate  vary  from  point  to  point  inside  a  tube,  the  relation 

between  T^and  jBy.  is  unique  for  laminar  flow  of  any  time 

D 

independent  fluid.  This  relationships  which  may  be 
conveniently  determined  in  a  cap i  1  lary= tube  viscometer 


defines  the  characteristic  behaviour 


a  fluid  just  as 


completely  as  shear  rate  -  shear  stress  relationship. 

In  general  n9  and  K9  may  be  expected  to  vary 


for  a  fluid.  However,  if  n5  and  K 5  are  constant 


i  .  e 


if  the  logarithmic  plot  of  T a n d  8v_  is  a  straight  line 

D 

=•  the  equation  (22)  becomes  the  power  law  for  a  fluid 
(equation  16).  Then 


n 


n 


n 


a  nd 


K  9  «*  K  <  — ) 


(22a) 


It  is  interesting  to  note  that  the  unique 

correlation  between  and  fy  is  applicable  only  for 

laminar  f  l  ovj  ,  Outside  the  laminar  flow  region  the  relation 

ship  between  X  and  is  dependent  on  the  pipe  diameter 

10  D 

(8,9) 

Equation  (22)  can  be  re-arranged  as  below  to 

define  the  friction  factor  f  for  laminar  flows 

A|p  D 


f  « 


L 


U3 


S3 
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sss 


\ 


64  k!  t'S)'0'-' 


<5  V1-"'  s"' 


whc  r  e 


64- 

r; 

kN  (8)n'-' 


(23) 


Comparing  equation  23  with  the  corresponding 
e q nation  for  N ewtonian  f 1 o w  vii 

£  m  ^  ^  -  —  -  —  -  -  ( e  q  u  n .  11) 

0 

Me  t  z  ner  (10,  1,2)  calls  K  a  generalised  Reynolds  nut.  be  r  . 

It  enables  the  data  for  both  Newtonian  and  non-Ne w t o n 1  a n 

fluids  to  be  represented  in  non-dimensional  form  in  a 

o 

plot  of  f  vs  R  .  For  laminar  flow  the  relation  between  f 

e 

o 

and  R  from  equation  (23)  would  be  a  straight  line  in  a 
logarithmic  plot*  Metzner  and  Dodge  (10)  have  extended 
the  idea  of  the  generalised  Reynolds  number  to  obtain  f 

9 

vs  Re  relation  for  turbulent  flow  also  as  in  Fig*  4* 


1.4.1  ANOMALOUS  FLOW  NEAR  A  WALL. 

In  the  derivation  of  the  equation  9  it  was 
assumed  that  there  was  no  "slip"  at  the  wa  1 1  *  However  in 
certain  non-Newtonian  fluids  (eg*  certain  high  polymer 
materials  and  certain  suspensions)  the  presence  of  a  boundary 
introduces  a  preferred  orientation  of  molecules  or  particles 
near  the  wall  in  an  otherwise  isotropic  material.  In  the 
case  of  suspension  of  solids  sometimes  a  thin  film  of 
lower  consistency  may  form  near  the  boundary  and  give  rise 
to  an  "effective  slip  velocity"  at  the  wall. 


17 
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Friction  factor  -  Reynolds  Number  Design  Chart 
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Oldroyd  (11)  calculates  the  "effective  slip 
velocity"  as  follows ; 

Assume  that  in  a  small  region  close  to  the  wall} 
i.e.  within  a  range  o<^<£  ,  the  shear  rate  will  differ 
from  (X)  » 


dt 


i-e-  ^  -  R-o 

Where  g  (T,D )  ~  0  when  ^  >  G. 


(24) 


The 


region  0<*y<£  is  called  the  region  of  anomalous 
The  velocity  u  just  outside  the  region  of 
anomalous  flow  is  then 

r  £ 

tw  -  K'O'b  ■+  J  1$  (5^)  ______ 

o 

;  wr  i  t  ing  j  =  S  (T) 


flow. 


(25) 


Ta.  -  s  (r) 


(26) 


S(T)  is  therefore  the  value  of  uat  yco  and  is  the  effective 
slip  velocity  at  the  wall.  It  depends  upon  the  local  shear 
stress  T * 


For  the  case  of  laminar  flow  through  a  round  tube9 


substituting  u(R)  =  s  (T^O 


<^v 

D 


4  s  CT*) 

R 


-V- 


in  the  derivation  of  equation  (9) 


Tu> 


t"  {tv  At 


(27) 


Where 


sv  

5  CT.)  - 
T  tTu)  = 


4-  "5  (Tu>) 

R 

S  (.  (uo) 

j  T^-fcO  A  T 


“  (28) 
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■5  CT“)  ^  called,  by  Oldroyd  (11),  the  effective 
slip  coefficient  and  is  the  effective  velocity  of  slip  per 
unit  shear  stress  at  the  vail. 

Equation  (28)  can  be  solveci,  in  principle,  for 
^(T)  and  f('T)  ,  by  making  observations  on  different  diameter 
tubes.  If  there  is  any  a noma lous  behaviour  the  plots  of 

C2W  _J _  rr 

T)  "p-  vs  will  be  dependent  on  the  value  of  D. 

However,  it  must  be  remembered,  that  -ppp  vs  lu*  curves 
will  be  dependent  on  D  at  high  values  of  X  due  to  onset  of 

bO 

turbulence.  This  should  not  be  misinterpreted  as  wall  effect. 


1. 4. 2  PARTICLE  SLIP  VELOCITY  IN  VERTICAL  PI P&S,. 

For  steady  flow  of  suspensions  in  vertical  pipes 
Kada  and  Henra  tty  (23)  observe  that  even  though  the  con¬ 
centration  at  a  section  be  constant  with  respect  to  time 
the  average  velocity  of  the  particles  could  be  different  from 
that  of  the  fluid.  The  difference  between  the  average  velocity 
of  the  particles  and  the  average  velocity  of  the  fluid 
is  called  " the  particle  slip  velocity”.  The  free  fall 
velocity  of  the  particles  is  a  close  approximation  to  the 
"particle  slip  velocity”  (23)  if  the  concentration  of  the 
suspension  is  not  appreciable. 

This  "particle  slip  velocity"  which  is  different 
from  the  "effective  slip  velocity"  of  Oldroyd,  is  important 
in  the  capillary  viscometer  studies.  If  the  volume 
concentration  of  solids  is  determined  from  the  quantity  of 
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suspension  discharged  from  the  capillary  tube  in  a  known 
interval  of  time,  the  observed  concentration  will  tend  to 
be  greater  than  the  actual  concentration* 


* 


&HAEIE&  XL 

FLOW  BEHAVIOUR  OF  SUSPENSIONS 


2*1  A  suspension  being  a  system  of  solid  particles 

dispersed  in  a  liquid  medium  one  can  visualise  many  types 
depending  on  the  physical  c ha r a c t e r i s t i c s  of  the  solid  and 
the  liquid  phases*  However*  only  those  suspensions  will  be 
considered  in  which  the  liquid  is  Newtonian  and  the  solids 
are  rigid,  do  not  react  with  the  liquid  and  do  not  possess 
electrical  charge* 

Generally  one  can  expect  the  "apparent  viscosity” 

v?  of  such  a  suspension  -  defined  as  the  ratio  of  shear 
va. 

stress  to  shear  rate  at  a  specified  shear  rote  -  to  depend 
on : 

(1)  the  concentration  of  solids 

(2)  the  shape,  size,  and  mass  of  the 

solid  particles 

(3)  the  viscosity  and  density  of  the 
suspending  medium 

( 4 )  the  s hear  r a  t  c 

Frisch  and  Simha  (12)  consider  the  suspension 
in  three  concentration  regions  wherein  a  characteristic 
behaviour  is  present,  viz, 

(a)  extremely  dilute  suspensions 

(b)  dilute  suspensions 

(c)  concentrated  suspensions 
Taking  these  in  detail: 

,2*1,  1  EXTREMELY  DILUTE  SUSPENSIONS 

In  these  the  suspension  is  so  dilute  that  the 
distances  between  the  particles  are  effectively  infinitely 
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great  in  comparison  to  their  sizes.  During  motion,  every 
particle  produces  a  disturbance  in  an  otherwise  laminar 
flow  around  it  in  a  restricted  field,  The  disturbances 
due  to  two  adjacent  particles  dc  not  overlap, 

Einstein,  in  1905,  found  a  mathematical  solution 
for  the  case  of  rigid  spheres  in  a  N e w t o  n i a n  l  i  q u i d ,  u n d  c  r 
t  h c  a  s  sump  t i on  o  ?  extreme  dilution  (  4 )  a  s 


who  r  e 


Yl  K  1  +2,5  C  (29) 

W 

’Y)  ~  r  e  1  a  t  i  v  e  v  i  s  c  o  s  i  t  y 

Vy 

-  4  p  pa  rent  v  i  :  c  o  s  1  t  y  o,f  ...lh&,  .s.u.S.llbn 
Viscosi  1  of  the  suspending  medium 

C  »  V  o  1  ume  c  otic  en  t  r  a  t  ion 


_  Vo  1  ume  o f  sol  i  d ,s 

V o 1 u to e  of  t  h e  s u s p e n s i o n 

Equation  29  indicates  that  the  suspension  behaves 
as  Newtonian  liquid  in  this  region  of  extreme  dilution,  as 
Yj  i s  i n d e pendent  of  the  shear  rate.  Einstein's  e q u a t i o n 

\y 

has  been  found  to  be  valid  for  values  of  concentrations  up 
to  about  3  percent  ( 4 ) * 

Following  Einstein’s  a  p  p  i*  o  a  c  h  ma  n  y  theoretic  a  1 
studies  have  been  conducted  on  infinitely  dilute  suspensions 
of  solid  particles  of  various  shape  (12)  and  the  result  can 

h e  genet a  1  1  y  e  x  p r c  s s e  d  as 


vhe  r  e 


1  -V-  a ,  C 


(  3  0 ) 


is  a  function  of  dimensionless  geometric 
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Daramcters  characterising  the  shape  of  the 
particle , 

It  is  interesting  to  note  that  the  densities  of 
the  solid  and  liquid  components  of  the  suspension  do  not 
have  any  effect  on  the  consistency  of  the  suspension  in 
this  range. 

Thus  the  infinitely  dilute  suspension  -  or 
real  suspension  with  volume  concentration  less  than  3  percent 
-  behaves  like  a  Newtonian  fluid  and  the  relative  viscosity  is 
n  linear  function  of  the  concentration* 

2^1 t2  DILUTE  SUSPENSIONS. 

In  these  suspensions  the  di sturbances  produced 
by  the  adjacent  particles  interact,  The  dependence  of 
viscosity  on  concentration  is  no  more  linear.  Numerous 
theoretical  and  empirical  solutions  have  been  proposed  by 
assuming  various  models  to  represent  the  interaction*  A 
good  review  of  works  in  this  range  is  given  by  Frisch  and 
Simha  (12),*  Generally  the  concentration  dependence  of 
viscosity  is  expressed  as  a  polynomial  with  the  values  of 
coefficients  depending  on  the  shape  of  the  particles* 

However  the  effect  of  shear  rate  is  not  considered  in 
any  of  the  work  reviewed  in  ref,  (12), 

2.1.3  CONCENTRATED  SUSPENSIONS,, 

In  these  the  particles  will  be  so  close  together 
that  they  form  networks,  with  flow  immobilising  power  far 
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greater  than  that  in  the  dilute  suspensions*  This  regime 
corresponds  to  a  concentration  in  which  particles  may 
form  close  packing,  Frisch  and  Simha  (12)  observe  that 
non-Newtonian  behaviour  is  to  be  expected  in  a  real 
suspension  in  this  range,  A  review  of  the  theoretical  works 
in  this  field  ;;iven  in  r e f  »  (  1 2 )  i n d  i c a t e  that  the  formulas 
derived  from  the  theoretical  approach  have  not  considered 
the  effect  of  shear  rate. 

Thus  the  theoretical  and  semi  theoretical  npproachc 
to  the  problem  of  viscosity  of  dilute  and  concentrated  sus¬ 
pensions  have  all  assumed  Newtonian  behaviour  and  other 
simplifying  pictures  of  particle  interaction,  The  result 
is  a  host  of  fomulne,  which  can  be  described  as  confusing, 
to  say  the  leant,  when  it  comes  to  the  application  to  a 
specific  case,  This  has  been  demonstrated  in  a  dramatic 
way  in  a  review  by  Rutgers  (16), 

2i1m1, 4  Rutgers  (16)  has  reviewed  2  80  references  in  the 

field  of  colloidal  physics  and  rheology,  and  96  formulas 
proposed  for  the  effect  of  concentration  on  the  relative 
viscosity  of  a  suspension  of  rigid  spheres.  All  these 
96  formulae  purport  to  give  ^  -C  relations  for  the 
Newtonian  range,  Other  than  Einstein's  formulas  (equation 
29),  which  is  valid  up  to  C “  0 , 0  ,  Rutgers  found  that: 
only  5  of  the  remaining  formulas  have  validity  over  the 
whole  concentration  range.  Also  only  two  of  these  five 
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formulae  are  easy  to  handle  and  are  based  on  some 
theoretical  argument*  They  are 
Mooney’s  equation: 


where 


*  C 

\  -  sc 


K 

S 


-2,76 
-  1,29 


K 


or 


s 


2,50 

1,40 


Vand  *  s 


t(y  "  relative  viscosity 
equation: 

4r  -  (1  -  A,C  +  A2C2)n 


(31) 


(32) 


where 


A1  »  1,10 
A2  -  0,97 
n  *  2,50 


or 


Aj  -  1,00 

A2  -  1,20 

n  -  2,50 


These  formulas  were  derived  for  suspensions  of 
rigid  spheres  in  Newtonian  liquid.  The  specific  gravity  of 
the  solids  was  the  same  as  that  of  the  liquids  Also  the 
suspensions  were  assumed  to  behave  as  Newtonian  fluids  - 
presumably  at  very  low  shear  rates. 

In  another  paper  Rutgers  (17)  has  reviewed  21 
investigations  on  the  suspension  of  rigid  spheres  in 
Newtonian  liquid.  The  spheres  were  of  the  same  specific 
gravity  as  the  liquid  and  of  sizes  0,3  to  400  microns* 

The  distribution  of  any  one  grain  size  was  such  as  to 
include  a  moderate  amount  of  grains  of  adjacent  sizes* 
Rutgers  (17)  observes  that  the  investigations  reviewed  by 
him  seem  to  indicate  that  the  suspensions  of  the  above 
type  behave  as  Newtonian  up  to  C  »  0*25* 
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Also  he  indicates  that  there  seem  to  be  several  volume 
concentrations  where  the  flow  behaviour  shows  changes 
e,g,  at  C  -  0,02  -  0,15  -  0,20  -  0,25  -  0,45  -  C  max. 

However,  a  study  of  the  data  used  by  Rutgers  (17) 
shows  that  most  of  the  experiments  were  conducted  either 
at  constant  values  of  fori  or  the  range  of  shear  rate 
used  was  very  small.  The  maximum  shear  rate  used  in  any 
experiment  was  1860  seel  With  this  in  view  much  weight 
cannot  be  attached  to  the  conclusion  that  Newtonian  behaviour 
occurs  up  to  C  ■  0,25  in  all  suspensions. 

From  these  reviews  (12,16,17)  one  can  conclude  that 
inspite  of  the  numerous  works  that  has  been  done  on  this 
topic  in  the  field  of  colloidal  science  and  rheology,  there 
is  no  full  appreciation  of  the  various  factors  affecting 
the  behaviour  of  suspensions.  This  is  especially  so  with 
regard  to  non-Newtonian  behaviour, 

Z*1  BE H A VXQ.UR._0 SilSPENSJO^ 

The  earliest  study  on  dilatant  behaviour  was  by 
Osborne  Reynolds,  His  original  explanation  given  in  1888 
is  probably  the  best  explanation  for  dilatancy,  (1)«  "He 
assumed  that  suspensions  of  solids  at  high  concentrations 
when  at  rest  consist  of  densely  packed  particles  in  which 
the  percent  of  voids  is  small  and  perhaps  minimum. 

Sufficient  liquid  is  present  only  to  fill  the  voids. 

The  motion  of  shear  of  such  fluids  at  low  shear  rates 
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requires  only  small  shearing  stresses  since  the  liquid 
lubricates  the  passage  of  one  particle  past  another*  At 
increasing  rates  of  shear  the  dense  packing  of  the  solids 
is  progressively  broken  up  and  since  there  was  only 
minimum  of  voids  the  liquid  present  is  insufficient  to 
fill  the  voids  and  lubricate  them*  This  causes  the  shear 
stress  to  increase  more  than  proportionately  with  the 
shear  rate*"  This  explanation  suggests  that  nil  sus¬ 
pensions  of  solids  in  liquids  should  exhibit  dilatant 
behaviour  at  high  solid  contentful), 

Freundlich  and  Roder  working  with  starch- 
water  and  quartz-water  systems  found  that  an  increase 
in  the  viscosity  of  the  system  sometimes  occurs  with 
increneps  in  the  she***  rate  (13),  This  was  also  termed 
dilatant  behaviour  as  it  was  similar  to  the  dilatancy 
described  by  Reynolds,  even  though  the  mechanism  may 
have  been  different.  Thus  the  term  "dilatant"  has  since 
been  used  for  all  fluids  which  exhibit  the  property  of 
increasing  apparent  viscosity  with  increasing  rates  of 
shear,  Metzner  and  Whitlock  (13)  have  found  that 
volumetric  dilation  (as  in  Reynolds  theory)  may  occur 
quite  separately  from  dilatancy  in  the  rheological  sense* 

In  the  Chemical  Engineering  field  the  non- 
Newtonian  behaviour  of  suspensions  has  long  been  recognized, 
Thomas  (18)  has  completely  investigated  the  behaviour  of 
Thorium  oxide  and  other  materials  of  particle  size  0*35  -  13 
microns  by  assuming  the  Bingham  plastic  law.  However, 
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regarding  the  suspensions  of  sand  size  particles  which 
exhibit  dilatant  behaviour,  there  is  no  similar  work  which 
has  tried  to  unify  the  effects  of  various  factors  - 
concentration,  particle  size  etc.-  on  the  consistency  in 
the  non-Newtonian  range* 

Various  investigators  have  tried  to  determine 
the  onset  of  dilatancy  in  suspensions  and  their  result  is 
summarized  in  the  following  Table  No*  1  taken  from  ref*  (13)* 


TABLE  1 


SUMMARY  OF  TESTS  ON  DILATANCY 


INVEST  I  GAT  01 

SOLIDS 

LIQUID 

VISCOSITY 

CENTIPOISE 

PARTICLE 

SIZE 

(MICRONS) 

DILATANCY 

OBSERVED 

AT  CONCN * 

MAX. 

SHEAR 

RATE 

SEC 

METZNER 

Ti02 

1,7,42 

0. 2-1,0 

30  -  47 

1000 

Glass 

beads 

18,42 

24  -  32 

possibly 

1300 

64 

42 

53  -  62 

it 

1300 

42 

85  -  105 

none 

700 

42 

24  -  105 

none 

1300 

FREUNDLIC1 

&  RODER 

Quartz 

1 

1-5 

42  -  45 

** 

Glass 

154,241 

3,4 

none 

639 

beads 

136,234 

4,10 

none 

639 

71,94 

10,30 

none 

639 

VAND 

Glass 

about 

beads 

80 

100-160 

none 

100 

. 
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The  Table  1  indicates  that  the  search  is  for  a 
unique  concentration  at  which  the  flow  changes  from 
Newtonian  to  dilatant  behaviour  irrespective  of  the  other 
properties  of  the  solid  and  liquid  phase*  It  is  the 
contention  of  the  author  that  there  is  no  apriori 
reason  to  assume  such  a  unique  concentration*  The 
beginning  of  dilatant  behaviour  must  be  functionally 
related  to  the  liquid,  solid,  and  the  flow  properties  of 
the  system  viz*,  density  and  viscosity  of  the  liquid, 
concentration,  size,  shape  and  density  of  the  particle 
and  the  shear  rate  and  its  variation  in  the  system* 

Only  by  dimensional  analysis  and  operation  of  dimension¬ 
less  terms  can  it  be  hoped  to  determine  the  transition 
zone  uniquely  -  in  terms  of  the  other  variables, 

1*3  BAGNOLD * S  APPROACH 

Bagnold  (19,20)  has  introduced  a  new  concept 
for  the  analysis  of  flow  of  solids  in  suspension*  Con¬ 
sidering  a  simplified  model  of  an  array  of  solid  spheres  he 
has  shown,  as  outlined  below,  that  as  a  result  of  momen¬ 
tum  transfer  between  successive  layers  of  the  suspension,  a 
repulsive  pressure  Ps  between  the  grains  and  a  shear  stress 
^  should  exist*  This  shear  stress  Ts  i s  additional  to  any 

residual  stress  T  due  to  momentum  transfer  within  the 

■f 

intergranular  fluid*  i*e*,  It  will  be  additional  to  that 
necessary  for  the  motion  of  the  fluid  only* 


.  .  i  ,  .  -wJ-'X';  e  •  r  ■ 
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Suppose  a  mass  of  rigid  elastic  spheres  with 
uniform  diameter  is  arranged  in  a  tetrahedral  -  rectangular 
pattern  which  gives  a  minimum  void  space  (corresponding  to 
volume  concentration  C*)*  Let  it  be  uniformly  dispersed  so 
that  distance  d  between  the  centres  is  increased  to  bd 
(Fig,  5), 


is  cu 


The  free  distance  between  two  adjacent  spheres 


Then  bd  »  a+d 
or  b  « 


+  1 


4-  l 


where 


n 

A  -  **  "Linear  concentration" 


(33) 


-  -  -  (Proof  in  App*  B) 


where  C*v  »  Maximum  possible  value  of  C  when  ail 
grains  are  in  static  contact. 


C*  ■  0*74  for  perfect  spheres 

■  0«65  for  natural  uniformly  rounded 
grains  (20) 

Bagnold  assumes  that: 

(a)  there  is  no  slip  between  the  solids  and  the 
1 i qui d 

(b)  a  steady  uniform  velocity  gradient 

prevails  throughout  the  mixture 


(c)  kinetic  energy  per  unit  volume  is  maintained 
constant  by  frictional  losses 

(d)  in  addition  to  average  velocity  u»  the  motion 
of  each  sphere  consists  of  three-dimensional 
oscillations  caused  by  succession  of  elastic 
collisions  between  spheres  in  faster  and  slower 
moving  layers 


' 
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Fig.  5  Definition  Sketch 


Fig.  6 


Correlation  between  G2  and  N  (Ref.  19) 
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Then,  according  to  Bagnold,  between  two  adjacent 

layers  A  and  B  shown  in  Fig*  5, 

/ 

SU  -  Kbd  ------  (64) 

dy 

where  K  is  a  constant,  Ccnsiclc-r  the  motion  of  grains  A 
relative  to  grains  B,  Let  the  number  of  collisions  made 
in  —  seconds  be  a  function  of  \  i ,  e «  ,  fn  (^\)  S\J  collisions 

bV  . T” 

occur  in  unit  time*  The  number  of  grains  in  unit  area  in 

each  layer  is  *  Each  grain  in  A  collides  at  some  angle 

(bdj5 

oc  with  the  grains  B,  thereby  experiencing  a  total  change 
of  momentum  2,-rn  Svj  CosoC  i n  the  y  direction  and  2.  rn  Su  Sin©c 
in  x  direction,  (m  “  mass  of  sphere)*  This  gives  rise  to 
a  repulsive  pressure 

fe,  =  TTTz  —  ^  2mStCosoC  -----  (33) 
b  d»  cl 

and  a  shear  stress 

Ts  -  lp&  oc  -  -  - (36) 

3 

Substituting  m  w  d  and  for  a  and 

in  equa t i on  35 , 

-  \<^  ^  “X  Fn  C.X)  d.  cos  ------  (37) 

and 

T,=  Sin  OC  - . (38) 

Where  is  a  constant* 

Bagnold  d i st ingui she s  two  regimes  of  f low •  Inertial 
(  Corresponding  to  high  rates  of  shear)  and  Viscous 
(Corresponding  to  low  shear  rates)*  The  theory  as  above  is 
for  the  inertial  range.  For  the  viscous  range  in  which 


33 


fluid  viscosity  predominates,  Bagno  Id  gives  a  similar 
theory  (19)  with  the  relation  between  shear  rate  and 
shear  stress  given  by: 


Where 


Tt  --  +  + {'(A)) 


clu 


(39) 


T  -  total  shear  stress  due  to  the  effect  of  fluid 

viscosity  as  modified  by  the  presence  of  grains 


Bagno  Id  conducted  experiments  in  a  rotary 
viscometer  with  spherical  grains  of  diameter  0#13  cm,  and 
density  the  same  as  that  of  the  suspending  medium,  The 
shear  stress  necessary  to  prevent  the  movement  of  the 
inner  cylinder  and  the  dispersive  pressure  on  the  inner 
wall  were  measured.  It  was  found  that  both'('art<3  ^  when 
expressed  as  dimensionless  groups 


Ts  ?s  cl 


T 

.2 


x  ip 
k,  d.’ 


*  A'f 

bear  a  single  valued  empirical  correlation  with  the 


dimensionless  group 


N 


o 

This  correlation  between  G4-  and  N  obtained  by 
Bagno  Id  is  shown  in  Fig*  6*  This  relation  gives. 


T 


o.oi3  ^  C—  ) 

s  ^  d.'y  ' 


for  inertial  region 


(40) 
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be  tween 


-  2,25 


d\j 


for  viscous  region  ------  (40a) 

The  transition  from  viscous  to  inertial  region  is 
G2  -  100  to  3000. 


It  will  be 


which  is  of  the  same 


noticed  that 

%  ?.  A* 


f  orm 


xv 

cl  .  C 

ns  • 


dL-C  >) 

V<? 

d  U* 


.  R 

e  * 

"  shear  Reynolds  number  used  in  problems 
connected  with  sediment  transportation 
in  open  channels* 


2*4  CRITIC  A  L  C  ONCE  NT R AT  1  ON , 

llll  M*l-‘  •«—*<"%» «  v.  **-  nWl.'Mi*  to*.  — '  •  •  -WT  -to  V*—."'-,—'-*' 

In  the  rheological  behaviour  of  suspensions  there 
seems  to  be  a  critical  concentration  at  which  there  is  a 
sudden  change  in  the  apparent  viscosity.  Frisch  and  Simha 
(12)  have  indicated  that  at  concentrations  near  close  packing 
one  might  expect  strong  and  abrupt  increase  in  the  apparent 
viscosity  of  the  suspension.  Schack  et  al  (14)  have  reported 
the  occurrence  of  a  critical  concentration  at  which  the  apparent 
viscosity  rapidly  rises  with  a  small  rise  in  the  concentration, 


. 
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This  critical  concentration  was  at  C  -  0.50  for 

quartz,  feldspar,  glass  and  calcite  of  sizes  5-50  microns. 

For  gypsum  Ccr  was  at  0 , 45  and  for  mien  it  was  0.35. 
Moreland’s  (15)  experiments  of  suspensions  of  coal 
(80-250  mesh  size)  in  mineral  oil  indicate  that  at  00.30 
the  apparent  viscosity  starts  increasing  at  a  very  rapid 
rate  with  the  concentration.  Rutgers  (7)  mentions  several 
critical  concentrations  at  which  the  flow  behaviour  may 
change,  It  is  interesting  to  note  that  "for  spheres  at 
00,25  rotation  of  the  spheres  may  begin  to  be  hindered  and 
00,45  interlocking  will  start"  (7). 

Thus,  it  seems,  depending  upon  the  shape  of  the 
particle,  there  exists  a  critical  concentration  at  which  an 
abrupt  change  occurs  in  the  rheological  behaviour.  Also  one 
may  expect  this  critical  concentration  to  have  values  between 
0.30  and  0,50* 
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EXfEfiirePTAi  INVESTIGATION. 

2*1  ;NTROp,Hn,tON 

A  special  investigation  of  the  rheological 
behaviour  of  sand-water  suspensions  was  undertaken  in 
the  Hydraulic  Laboratories  of  the  University  of  Alberta* 

The  justifications  for  the  investigations  were: 

(1)  The  inadequacy  of  previous  work  (sec*  2.1  and  2*2) 

(2)  The  need  to  understand  the  simpler  aspects  of  concentrated 
sand  transport  in  water  as  a  preliminary  to  explaining  the 
nature  of  transport  in  the  saltating  cloud  in  the  duned 
sand  bed  of J a  canal  or  a  river*  This  later  unsolved 
problem  has  assumed  new  Importance  in  the  International 
Hydrologic  Decade* 

(3)  Possible  relation  of  the  results  to  means  for  estimating 
transport  in  pipes  from  flow  parameters* 

The  aims  of  the  investigations  were: 

(a)  To  develop  a  viscometer  for  the  rapid  and  yet 
accurate  determination  of  the  rheological 

i 

properties  of  the  suspensions* 

(b)  To  obtain  data  on  sand-water  suspensions  from  the 
above  viscometer  and  analyze  them  in  the  light  of 
present  knowledge  of  the  flow  of  non-Newtonian 
fluids  and  the  behaviour  of  suspensions* 

(c)  To  study  the  applicability  of  Bagnold’s  hypothesis 
to  laminar  flow  through  vertical  pipes* 
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3  ■.  2  CAPILLARY  VISCOMETER  , 

Many  instruments  considered  as  multi-speed  or 
multi  pressure  viscometers  are  available  commercially,  but 
most  are  rotating  concentric  cylinders,  rotating  discs,  or 
cone  and  plate  type.  They  are  usually  restricted  in  designed 
working  range  which  cannot  be  altered  easily.  Also  a  good 
rotating  viscometer  to  cover  a  wide  range  of  shear  values  with 
reasonable  accuracy  is  quite  expensive*  From  the  consideration 
of  these  points  it  was  decided  that  a  capillary  viscometer  be 
designed  and  built  to  cater  to  the  special  needs  of  testing  of 
sand  suspensions*  This  has  a  special  advantage  that  the  in¬ 
formation  obtained  from  such  an  instrument  can  be  directly 
related  to  flow  through  pipes  under  appropriate  conditions* 

The  essential  features  of  a  capillary  viscometer  are** 

(a)  A  Storage  Chamber*  This  should  have  an  agitation 
mechanism  for  preparation  ©f  suspensions. 


(b)  Device  for  precise  control 
pressure  * 

(c)  A  Capillary  tube* 

(d)  Discharge  measuring  device 

(e)  Temperature  Control, 


measurement  of 


There  are  different  capillary  viscometers,  with 
different  modifications  to  suit  the  particular  needs  of  the 
problem,  mentioned  in  the  literature  on  viscometry  (1,  21,  22, 
14) «  Noting  the  advantages  and  disadvantages  of  many  of  these, 
and  with  the  object  of  developing  a  simple,  convenient,  rapid 
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and  yet  sufficiently  accurate  instrument  the  following  final 
design  of  a  capillary  viscometer  was  built  after  many  trials* 
Fig*  7,  8,  shows  the  schematic  diagram  of  the 
viscometer  setup*  Figures  9  to  13  show  the  various  details* 

3*2*1  STORAGE  CHAMBER 

A  Chamber  of  capacity  500  c»c*  was  made  by  using 
a  standard  plastic  pipe  3M  O.D,  5  5/8"  long,  and  connecting 
the  two  ends  to  brass  plates  by  water  tight  seals  as  shown 
in  Fig*  9*  A  concentric  3%"  O.D*  plastic  pipe  provided  the 
water  jacket*  An  agitator  mechanism  consisting  of  two  sets 
of  blades  mounted  on  a  pair  of  brass  shafts  was  powered  by 
a  1/30  H.P*,  1,550  r.p.m*  motor  through  a  chain  drive. 

3*2*2  PRESSURE  REGULATIONS. 

The  mixing  chamber  was  connected  to  the 
pressurized  air  supply  through  a  0-30  p*s+i»  fine  control 
pressure  regulator.  The  pressure  in  the  chamber  was  measured 
by  a  manometer*  Bromoform  (sp*  gr *  2*87)  or  mercury  was 
used  as  manometer  liquid  depending  of  the  order  of  pressure 
to  be  measured* 

3*2*3  CAPILLARY  TUBE . 

The  capillary  tube  for  use  in  this  instrument  was 
formed  as  a  separate  attachment*  Two  sizes  of  capillaries 
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Fig.  7  Schematic  Diagram  of  Capillary  Viscometer 


Fig.  8  Schematic  Diagram  of  Constant  Temperature 
Recirculation  System 
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[0 

0 

0 

0 

-Flonged  Sealed 
Roller  Bearing 


1:1  Chain  Drive  from 
1/30  H.P. ,  1550 
RPM  Motor 


(0 

0 

0 

Rubber  Sealing  Washer- 
(cemented  to  shaft) 


1/4"  Brass  Plate 


3/8"  Plastic 
O-  Ring  Seal 

Wall  of  Inner  Vessel 


Fig.  10  Schematic  Diagram  of  Bearings,  Seal  &  Drive 


Standard  Poly  Flow 
3/8"  Brass  Fitting 


Tygon  Coupling 


Stop  Cock 


Water  Outlet  Connected  to 
Inlet  of  Mixer 

Water  Jacket 

Capillary  Tube 


Fig.  11  Details  of  Capillary  Tube  Connection 
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Fig.  13  THE  STORAGE  CHAMBER 


Fig.  12  THE  EXPERIMENTAL  SETUP 

(A)  STORAGE  CHAMBER 

(B)  CAPILLARY  TUBE 

(C)  MANOMETER 

(D)  PRESSURE  REGULATOR 


Fig.  16  PHOTOMICROGAPHS 

A  60  -  70  MESH  SAND 
B  70  -  80  MESH  SAND 
C  SCALE  IN  MM. 
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were  used  in  the  present  investigation.* 

Tube  Aj  122  cm,  long  2,0  ttttv  dia, 

B:  118  cm,  long  2,7  mm,  dia. 

The  tubes  were  of  glass.  Each  tube  had  a  stopcock 
fused  at  one  end  and  the  entire  length  of  tube  was  encased  in 
a  plastic  water  jacket  and  was  mounted  on  a  frame.  The 
capillary  tube  was  connected  to  the  bottom  of  the  storage 
chamber  as  shown  in  Fig,  11*  The  tube  was  set  in  perfect 
vertical  position.  The  connection  to  the  chamber  and  the  tube 
was  such  that  the  tubes  could  be  changed  easily  without 
disturbing  the  general  setup, 

3,2,4  DISCHARGE  MEASUREME NT, 

The  discharge  of  the  suspension  through  the 
capillary  at  the  lower  end  of  the  tube  was  at  atmospheric 
pressure.  The  flow  was  allowed  t©  establish  and  then 
about  50  c,c,  of  the  discharge  was  collected  in  a  known 
interval  of  time  in  a  container#  The  time  was  measured  with 
an  electrical  stop-watch  of  accuracy  1/10  sec* 

hJUl  TEMERAIiJBE  C0£im 

A  constant  temperature  recirculating  system,  as  in 
Fig,  9,  was  used  to  maintain  the  temperature  of  the  liquid  in 
the  chamber  at  desired  level.  A  thermostat  controlled  heating 
coil  maintained  the  constant  temperature  in  the  water  tank 
and  a  small  pump  circulated  this  water  through  the  capillary 
tube  jacket  and  the  mixing  chamber  jacket.  This  arrange¬ 
ment  gave  a  control  of  the  temperature  of  the  liquid  in  the 
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chamber  to  an  accuracy  of ±  1°F  • 

3*2*6  OPERATION* 

The  operation  of  the  capillary  viscometer  is 

relatively  simple*  Once  temperature  is  well  under  control 

the  only  readings  to  be  taken  are  manometer  reading,  the 

average  static  head,  and  the  volume  of  liquid  issued  in  a 

known  period  of  time*  If  the  capillary  length  and  the 

diameter  are  accurately  known  the  elements  of  the  flow 

curve  <T~  _  D  and  & v 

4L"  D 

can  be  calculated* 

However,  it  is  very  difficult  to  measure  the 
diameter  of  a  capillary  accurately*  Also,  in  the  dilatant 
non-Newt onian  fluid  the  pressure  difference  in  a  capillary 
tube  is  very  sensitive  to  the  diameter*  Thus  there  is 
always  a  possibility  of  an  error  in  the  flow  curve  due  to 
an  error  in  the  measurement  of  the  diameter* 

Also,  the  measured  pressure  difference  will  have 
to  be  corrected  for  (a)  loss  at  entrance  (b)  loss  at 
the  valve  and  (c)  the  kinetic  energy  at  the  exits  But 
the  magnitudes  of  these  losses,  especially  in  a  non~ 

Newtonian  fluid  flow,  are  difficult  to  predetermine  and 
corr ec  t  * 

To  overcome  these  difficulties  an  accepted 
method  is  to  calibrate  the  instrument  with  stable  Newtonian 
liquids  of  known  viscosity*  The  calibration  will  be  more 
meaningful  if  the  ranges  of  shear  rate  and  shear  stresses  used 
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in  the  calibration  are  comparable  with  those  used  in  the 
investigation*  The  calibration,  as  discussed  below, 
achieved  this  purpose  quite  satisfactorily* 

1*2*1  QAIMMXIQILl 

Two  oils,  M  and  N,  whose  exact  viscosities  at 
various  defined  temperatures  were  known  accurately,  were 
used  to  calibrate  the  two  capillary  tubes  A  and  B}  The. 
viscosities  of  the  oils  were  as  follows; 


TABLE;  2 

VISCOSITIES  OF  THE  CALIBRATING  OILS 


VISCOSITY 

IN  C.P.  AT 

70°F 

90°F 

100°F 

Oil  M 

4.532 

3.372 

2.939 

Oil  N 

9.162 

6.274 

5.304 

The  flow  rates  at  various  pressure**  were 
measured  (Tables  A- l  to  A~4)  and  the  flow  curve  plotted 


as  in  Fig*  14,  and  Fig*  1.5*  In  these  figures  X « 


torn 


4  L 


viher  e 


"■  wall  shear  stress  computed  from 


A|p  «  observed  pressure  difference 
L  *  measured  length  of  the  capillary 
D  •  measured  diameter  of  the  capillary 


From  the  plots  in  Fig*  14  and  15,  it  is  found  that 
the  observed  viscosity  is  greater  than  the.  actual  viscosity* 


. 
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Fig.  14  Calibration  of  2  mm  Tube 
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By  considering  a  correction  to  be  applied 


to  observed 


shear  stress  only 


T 


tom 


cC 


where 


corrected  vail  shear 


stress 


(41) 


and  oC  **  a  coefficient  to  take  care  of  various 

pressure  losses  and  other  inaccuracies 
in  the  measurement « 


It  is  found  that  oC  i s  a  constant  for  a  given  tub® 

for  the  range  of  viscosities  of  oil  M  and  oil  N* 

Hence,  a  value  of 

oC  »  0*82  for  tube  A 
and  oC  «  0*90  for  tube  B 

is  adopted  as  a  result  of  the  calibration*  These  coefficients 
are  used  in  the  calculation  of  shear  stress  of  the  suspensions 
tested  in  the  capillary  viscometer. 


2*2*8  ISL  2M&  suspensions^ 

The  sand  used  in  the  tests  was  Ottawa  sand,  passing 
through  #60  sieve  and  retained  on  #80  sieve*  Fig*  16  shows 
the  photo-micrographs  of  the  sand*  The  average  size  of  d" 
0*22  mm  was  accepted  as  representative  and  was  used  in  all 
the  calculations. 

The  object  of  the  tests  was  to  obtain  the  flow 

SV 

curve  elements  (  <T  and  u  )  for  various  concentrations, 

u)  U 

over  as  wide  a  range  of  shear  stress  as  possible  with  the 


The  tests  covered  a  range  of  concentration 


instrumen  t « 


■ 
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from  C  661  0 •  1  .3  to  C  **  Q#35#  All  tests  were  performed  at  a 
constant  temperature  of  70°F* 

In  a  typical  test,  corresponding  to  a  concentration 
Clt  the  quantity  of  dry  sand  and  water  required  to  make 
450  c«e®  of  suspension  of  concentration  was  determined 
gravimetr ically «  The  specific  gravity  of  sand  was  assumed  as 
2*64*  The  required  quantity  ©f  dry  sand  was  placed  in  the 
filler  and  flushed  in  to  the  chamber  with  the  previously 
calculated  quantity  of  water*  The  mixture  was  brought  to  a 
suspension  by  turning  on  the  agitator*  The  temperature  was 
checked  by  drawing  out  a  small  quantity  of  the  suspension  by 
the  drain  valve#  When  the  required  constant  temperature  of 
70°F  was  established,  the  chamber  was  pressurized  to  give 
the  required  shear  stress#  The  capillary  valve  was  then 
opened  and  some  suspension  allowed  to  drain-out  to  establish 
the  t  low*  Then  about  50  c»c*  of  the  suspension  was  collected,* 
The  exact  volume  of  suspension  collected  and  the  time  taken 
were  noted#  Also  the  weight  of  the  suspension  was  noted  and 
the  suspension  dried  in  an  oven  and  the  weight  of  dry  sand 
left  was  noted*  From  these  measurements  the  exact  concent¬ 
ration  and  the  volume  rate  of  flow  were  calculated#  It  was 
found  that  in  the  instrument  the  concentration  intended  and 
the  exact  concentration  obtained  in  a  test  differed  within 
■Jl0*Q2„  This  was  not  considered  as  a  serious  defect  as  the 
object  was  to  obtain  the  flow  curve  data  over  a  range  of 
concentrat ions* 
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Tables  A-5  and  A-6  give  a  summary  of  the  data* 

From  the  observed  values  of  pressure  drop  the  wall 
shear  stress  T  is  calculated  as  T  s  oC  X  .  *  The  flew 

U)  viJo  com 

function  8  V  is  calculated  from  the  discharge  measurement? 

D 

Tables  A  -  7  and  A  -  S  list  the  values  of  “T  and  8  V  and  other 

calculated  quantities  corresponding  to  various  experiments? 


CHAPTER  IV 


INTERPRETATION  OF  DATA 

4.1  &HBM  il&m  *  FLOW  FUNCTION  RELATIONSHIP 

Flov  curves,  i*e*,  Plots  of  wall  shear  stress 

against  flow  function  defined  by  8V  ,  are  the  basis  for 

D 

studying  rheological  behaviour  from  capillary  viscometer 

data*  Fig*  17  shows  a  plot  of  shear  stress  and  the 

flow  function  S2L  »  on  an  arithmetic  scale*  The  concentration 

D 

C  Is  used  as  the  third  parameter*  All  the  data  obtained  from 

both  the  tubes  A  and  B  are  plotted*  By  Interpolation  among 

plotted  points,  lines  representing  the  relationship  between 

T  and  ££  ,  and  having  a  constant  value  of  C,  are  drawn* 
w  D 

Five  such  lines  having  a  value  of  C  ■  0*13,  0*19,  0*26,  0*29 
and  0*34  are  shown  in  Fig*  17*  It  can  be  seen  that  these 
lines  are  curves  with  concave  curvature  towards  the  'TL, 
axis  and  thus  indicate  the  non-Newtonian  behaviour  of  the 
suspension* 

As  a  first  approximation,  it  was  assumed  that 
at  very  low  shear  rates  the  fluid  behaviour  can  be  approximated 
to  that  of  a  Newtonian  fluid*  With  this  assumption  the 
relative  viscosities  of  the  suspension,  corresponding  to 
various  concentration  lines,  were  calculated  at  the  lowest 
shear  rate  up  to  which  the  lines  can  be  reasonably  ex¬ 
trapolated*  These  values  were  then  compared  with  the 
formulas  of  Vand  (equation  32)  and  Mooney  (equation  31) 
which  have  been  recommended  by  Rutgers  (16)  as  applicable 
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Fig.  17  Plot  of  rw  v*  6V/D  (Arithmetic  Scale) 
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to  the  entire  concentration  range  of  the  suspension  of 
rigid  spheres  behaving  as  Newtonian  fluid  (sec*  2*1)* 

Even  though  Einstein's  (equation  29)  is  not  applicable 
to  the  present  case,  it  was  used  to  calculate  the  relative 
viscosities  to  get  an  idea  of  the  variation  of  values 
obtained  by  the  other  two  equations  and  from  the  Fig*  17* 

Table  3  gives  a  summary  of  the  relative  (to  water)  viscosities 
calculated  as  above* 


RELATIVE  VISCOSITIES  OF  SUSPENSION 


Concn* 

Relative 

Viscosities  by  Formula 

Relative  Viscosity 

C 

VAND 

MOONEY 

EINSTEIN 

crom 

tig.  17 

0*13 

1*328 

1*480 

1.325 

• 

1.85 

At  Si  m 

0.19 

1*480 

1.920 

1.475 

2.50 

D 

0.26 

1.645 

2.770 

1.650 

3.30 

2000  Sec."l 

0.29 

1*780 

3*380 

1.750 

4.25  . 

0*34 

2*080 

5*380 

1.850 

16.70  - 

-At  £L  m 

D 

1000  Sec."*1 

^  of  water  is  assumed  as  1 

centipoise 

Even  though  this  comparison  is  very  approximate 
it  clearly  Indicates  that  the  above  formulas  are  not 
applicable  to  the  conditions  of  the  present  investigation* 

The  pronounced  effect  of  concentration  on  the  characteristics 
of  the  suspension  is  apparent  from  the  Fig*  17* 
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Since  the  behaviour  of  the  suspension  is  non- 

Newtonian,  a  better  method  of  representing  the  Capillary 

Viscometer  data  is  to  plot  *£*  vs  8V  in  a  logarithmic  plot* 

W  D 

It  can  then  be  analyzed  in  the  light  of  the  non-Newtonian 
fluid  behaviour  discussed  in  sec*  1*3* 

iLt2.il  ANOMALOUS  FLOW 

From  the  nature  of  the  data  collected  in  the 
present  investigat ion,  it  was  not  possible  to  apply 
Oldroyd's  method  (sec*  1*4*1)  to  determine  effective 
slip  velocity  at  the  wall*  It  is  assumed  that  there  is 
no  anomalous  behaviour  of  the  suspension* 

4.2.2  PARTICLE  SLIP  VELOCITY 

The  particle  slip  velocity  (sec*  1*4*2)  can 
approximated  to  the  free  fall  velocity  at  low  concentration 
of  solids*  The  free  fall  velocityVIf  of  the  sand  size 
d  ■  0*22  mm  is  2*80  cm/sec*  In  the  concentration  range 
C  ®  0*13  to  0*30  the  velocities  of  the  suspension  in  the 
tests  were  between  60  cm/sec*  to  210  cm/sec*  Thus*  if 
it  i 8  assumed  that  the  particle  slip  velocity  in  this 
concentration  range  is  2*80  cm/sec**  it  will  be  a  small 
fraction  of  the  velocity  of  the  suspension*  Hence*  it  is 
assumed  to  be  too  small  to  affect  concentration  measure¬ 
ments* 

However*  in  the  concentration  range  C  -  0*30  to  0*35 
the  velocities  of  the  suspension  were  between  10  cm/sec*  and 
40  cm/sec*  From  Fig*  17  it  can  be  seen  that  there  is  a 
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pronounced  effect  of  concentration  of  solids  which  has 
Increased  the  relative  viscosity  of  the  suspension  to 
as  high  a  value  as  16*  It  is  therefore  felt  proper  to 
assume  the  particle  slip  velocity,  at  this  high  con¬ 
centration  range,  to  be  given  by  the  hindered  settling 
velocity.  The  following  empirical  formula  has  been 
recommended  by  Badger  and  Banchero  (24)  for  the  calculation 

of  hindered  settling  velocity  of  spherical  particles  of 

-  / 

uniform  size; 


wh 

where 

Wf 

£ 


Wf  'I'CO  . (42) 

hindered  settling  velocity 

free  settling  velocity 

volume  concentration  of  voids  -  1  -  C 


*  -  1.82  (1-6) 

CO  •  £2  X  10  for 

values  of  £between  0.95  and  0.50* 

Assuming  this  formula  to  be  applicable  to  the 

sand  used  in  the  tests,  the  hindered  settling  velocity 

of  the  sand  at  C  ■  0*33  is  calculated  as; 

For  C  -  0.33 

6  -  0.67 

2  -  1.82  (0.33) 

^<0-  (0.67)*  x  10.  _  0#ll3 


W 


2.80 


Wh  -  2.80  x  0.113  -  0.316  cm/sec. 

Considering  this  as  the  representative  particle 
slip  velocity  in  the  concentration  range  C  -  0.30  to  0.35, 
it  will  be  noticed  that  the  particle  slip  velocity  is 


\ 


7 
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between  3%  to  0*757*  of  the  velocity  of  the  suspension  in 

this  range*  Hence  it  is  assumed  too  small  to  affect  the 

measured  concentration  values* 

AjJL  LOGARITHMIC  FLOW  CURVE 

Fig*  18  is  a  logarithmic  plot  of  against  ££ 

D* 

All  the  data  shown  in  Fig*  17  are  plotted  but  the 

concentration  values  have  been  stated  at  intervals  of  0*03 

and  then  marked  distinctively*  It  can  be  seen  that  for  a 

given  value  of  concentration  the  variation  of  ^with  ££ 

D 

can  be  represented  by  a  straight  line*  Even  though 
there  is  a  slight  scatter  of  data  the  trend  is  very  clear* 
For  concentrations  up  to  0*30  the  lines  of 
constant  C  are  parallel*  Thus  they  can  be  expressed  as 


The  value  of  K'  depends  upon  the  value  of  C*  It  is 
found  that  n*  -  4/3,  indicating  that  the  behaviour  of 
the  suspension  is  dilatant  in  the  range  of  shear  rates 
tested* 

In  the  concentration  range  C  -  0*30  to  0*35  * 
there  is  a  sudden  change  in  the  characteristics  of  the  fluid* 
The  line  C  -  0*35  is  still  straight  but  the  value  of  n9  is 
1*04*  The  value  of  K9  has  undergone  a  sudden  Increase* 

The  lowest  concentration  tested  was  0*13*  The 
non-Newtonian  behaviour  is  evident  at  this  concentration 
also  for  shear  stress  values  greater  than  50  dynes/cm^ 
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Fig.  18  Plot  of  rw  vs.  8V/D  (Logarithmic  Scale) 
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and  flow  function  value  greater  than  2500  sec-1.  Thus 
the  conclusion  of  Rutgers  (17)  that  25%  volume  concent¬ 
ration  as  the  minimum  for  the  onset  of  non-Newtonian 
behaviour  of  suspension  of  rigid  spheres  is  not  true 
for  the  sand  suspensions  tested. 


For  values  of  C  between  0.13  and  0.30  the 


shear  rate 


The  minimum  and  maximum  shear  rates  tested  in  this 


-(Equn.  21) 


.  -I 

and  G500  sec  • 

concentration  region  are  therefore  2300  sec-1^  There  is 
definite  and  unmistakable  proof  of  dllatant  behaviour  at 
all  concentrations  between  0.13  and  0.30  in  this  range  of 
shear  rates.  This*  when  compared  with  the  maximum  shear 
rate  of  1300  sec"1  and  the  conclusions  of  Table:  1,  show 
the  importance  of  shear  rate  value  on  the  onset  of  non- 
Newtonian  behaviour  of  suspensions. 

Fig.  18,  even  though  useful  in  indicating  the 
non-Newtonian  behaviour  and  the  effect  of  concentration, 
is  not  a  dimensionless  plot.  So  is  unlikely  to  give 
simply  the  possible  relative  influence  of  various  parameters 
of  a  suspension.  To  obtain  potentially  meaningful  dimensionless 
terms  Bagnold's  hypothesis  was  adopted  as  a  guide. 


JuA  ADAPTATION  OF  BAGNOLD'S  CONCEPT 

The  essence  of  Bagnold's  concept  (Sec.  2.3)  is 
the  recognition  of  the  shear  and  direct  stresses  due  to 
impact  of  grains  and  its  unique  correlation  with 
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?  and  d*  However,  his  experiments  were  conducted  with 

uniform  shear  strain  throughout  the  flow,  l*e«,  ,£u  *  constant, 

dy 

and  with  specific  gravity  of  the  spheres  same  as  that  of  the 

fluid*  In  the  present  investigation  the  velocity  gradient  is 

constant  across  the  tube,  the  particles  are  not  spherical 

and  -  2*64  gm/cm^  and  -  1*00  gm/cm^* 

In  the  following  analysis  the  concentration 

range  C  -  0*13  to  0*30  is  considered  as  the  range  of  interest* 

The  same  analysis  is  however  applied  to  the  range  C  ■  0*30  to 

0*35  even  though  the  behaviour  in  this  range  may  be  different* 

In  Fig*  18,  it  has  been  shown  that  Twand  JgX.  follow 

D 

the  equation 


T  -  K»  (•22. )n 


(22) 


and  that  n*  is  constant  at  4/3  for  various  concentrations  up 
to  0*30*  Hence,  in  the  range  of  shear  rates  tested  the 
•power  law*  and  equation  16  are  applicable* 
n  *  ■  n  -  4/3 


k  '  -  k(<2b±1) 
4n 


n 


(22a) 


K  -  fn  (C) 


From  equation  17,  v  /  4** 


and  hence  the  average  shear  rate  is  proportional  to 

D 

Also  the  maximum  shear  rate  is  given  by 


v 


i  «  e  • 


max 


"dm 


) 


w 


3p+l 

4n 


(-  — ) 
v  ^  L 


o<^ 


8V 

D 


'  D  1 


not 
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Thus  the  average  shear  rate  and  the  maximum  shear  rate 
are  proportional  to  §£  *  It  follows  therefore  that  the 


term  adequately  represents  the  shear  rate  in  the  tube* 
D 

Suppose  the  total  shear  stress  in  the  capillary 


is  considered  to  be  made  up  of  two  parts 

T*  -T£  +T 

r  s 


(44) 


Where  ^  s  -  Shear  stress  due  to  the  presence  of  solids  in 
the  suspension  as  per  Bagnold*s  concept* 

-  Shear  stress  due  to  the  flow  of  the  liquid 
medium  alone* 


As  a  first  approximation  assume 

Tf  -  TJ(tt) 


(45) 


Where  m  -  coefficient  of  dynamic  viscosity  of  the 
Newtonian  liquid 

Ts  -Tw  -v)  (i2L) 

L  D  -  -  -  - 


(45a) 


Now  the  grain  shear  stress  should  depend  upon  the 


characteristics 

of 

(a) 

Solids  in  suspension:  sice  d ,  shape* 
f g*  concentration  C* 

density 

(b) 

Liquid 

suspending  mediums  denslty?f» 

vi scosi 

(c) 

Flows 

the  flow  function 

D 


Ts  “Fn  (d,Pg,  C,  Shape  factor,  Pf,f]  , 

-  - - (46) 

The  linear  concentration  X  is  considered  as  a  better 
representative  parameter  of  the  crowding  effect  of  the  particles 
than  the  volume  concentration  C*  Replacing  C  by  X  » 


as  has 


« 


■ 


It  a  up 
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been  done  by  Bagnold, 

Ts  *  Fn(  j^L,  Shape  factor) 

s  1  D 

- - (47) 


By  dimensional  analysis. 


Ts  f ,  df 


-  fn 


Shape  factor 


(48) 


These  dimensionless  groups  must  be  functionally  related,  if  the 
hypothesis  is  correct*  In  the  present  investigation  fyyfs  and 

<r  P 

shape  factor  were  constant  for  all  experiments*  So,  - 5 — 


with  X  as  a  third  parameter. 


is  plotted  against (.g3£.\  ^  ^ 

D  1 

in  Fig.  19.  it  can  be  seen  that  for  a  given  value  of  X  , 


there  is  a  good  correlation  of  the  two  parameters*  The 
relationship  is  linear  for  most  part  of  the  plot  l*e*  for 
values  of  (  ^  up  to  600* 


Thug  for  a  constant  value  of  X  In  this  region 


%  %  dl 

T 


2.  .» 


*1  (  83L:is  ci  ^ 
v  D  r[ 


(49) 


The  value  of  m  ■  5/3  for  X  <  3*5 
and  m  ■  1*04  for  X  >  4*0 


Also  Kj  -  fn  ( X ) 

Table:  4  shows  the  value  of  calculated  from  Fig*  19 
for  various  value  of  X  *  To  find  the  functional  relation 
between  Kj  and  X  ,  Kj  is  plotted  against  X  in  Fig.  20.  For 
values  of  X<  3*5,  there  is  a  clear  indication  that 
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(8V/D)/>,d% 

Fig.  19  Variation  of  rt/>td2/i^  with  (8V/D)/9$d  /tj  and  X 
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TABLE  4 


VALUES  OF  Kj  FOR  VARIOUS  VALUES  OF  X 


•<1 

2.35*  102 

2.85  x  10‘2 

3.25  *102 

4.65x10  2 

X 

1.4 

1.8 

2.0 

3.0 

Fig.  20  Variation  of  Kj  with  X 


0.06 


■  : 
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*1  -  K2  .  X 


(50) 


Where  K2  is  a  constant*  From  fig*  20  K2  •  0*016* 

Thus  for  \<  3*5*  and  ^  and  shape  factor  being 


constant 


T.  A1 


-  K. 


av 

D  ^ 


,z  -  ® 


] 


(51) 


It  is  interesting  to  note  that 


ts  ?s  a2 

m1 


is  the  sane  as  the  parameter  Q  used  by  Bagnold*  From 

T  2  z 

equation  51  it  follows  that  a  plot  of  G)  against  8V  d 

T  D 

must  give  a  unique  correlation  for  all  values  of  X 


less  than  3*5* 


Fig*  21  is  such  a  plot*  All  the  experimental 
points  are  represented*  It  will  be  seen  that  there  is  a 
good  correlation  of  the  data*  All  the  data  with  X<3.5 
fall  on  a  single  curve  -  (Curve  A)*  Even  though  a 
straight  line  as  per  equation  5 1  can  be  drawn  through  these 
polnts&  a  smooth  curve  as  In  Fig*  21  fits  the  data  very  well* 
The  data  with  X>  3*5  are  in  a  different  regime* 

The  points  with  vqlues  ©f\>4*0  seem  to  follow  a  relation 
similar  to  those  with  X<3*5* 

It  is  interesting  to  observe  that  the  values  of 

2 

G,  of  all  the  data  fall  between  180  to  1500o  This 
T 

corresponds  to  the  transition  region  between  viscous  and 
inertial  zones  of  grain  shear  indicated  in  fig*  6o 

Consider  the  curve  A  (X<3«5)*  Even  though  the 


2,000 


65 


-  o  o 

o  o 

o 


zux/z^ 


66 


abscissa  is  different  between  this  and  curve  in  Fig*  6 ,  the 

similarity  is  striking*  Both  the  curves  depict  the  variation 

of  grain  shear  stress  with  the  shear  rate,  but  the  nature  of 

the  shear  is  different*  For  the  same  range  of  the  ordinate 
2 

the  shape  of  the  curve  is  similar*  Considering  both 

the  curves  as  similar,  curve  A  can  be  considered  as  a 

transition  between  a  line  of  slope  1*^0  and  another  of 

2 

slope  2*00  as  shown  in  Fig*  22*  At  G  the  curve  changes 

T  cr 

to  a  straight  line  of  slope  1*0* 

For  value*  of  cr 

x^ 

t,  =  4M  iV) 

Tw  =  T,  +  Ts  =  fi  +  k*X  ] 

(  L  J  [  D  -  - - (52) 

which  Indicates  a  Newtonian  flow  behaviour*  Similarly  for 

all  ^  y  cr  9  since  the  slope  of  the  line  in  Fig*  22 

will  be  greater  than  unity  the  rheological  behaviour  would 

be  non-Newtonian  dilatant* 

Thus  by  considering  the  curve  A  to  be  similar  to 

Fig*  6,  a  critical  value  of  the  term  i*e*  G||  can  be 

envisaged  such  that  the  rheological  behaviour  is  Newtonian 

2  2 

for  values  of  G,  ^  For  values  of  X<3*5  this 

1  i  C7 « 

cr  may  be  expected  to  depend  upon  ^ /?s  also  the 

shape  factor  of  the  solids  in  suspension.  For  \>  3*5  a 


21 
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similar  relation  may  exist,  bur  the  experimental  data  is  scanty 


to  warrant  any  conclusions* 


Since  the  relation  between  G|T  and 


U’ 


D  Y! 


is  unique  it  follows  that  the  term 


D 


will  also 


2  -  ^ 

be  critical  at  <3)T  cr*  This  indicates  that  the  change  from 

Newtonian  to  non»Newtonian  behaviour  occurs  at  a  unique 

combination  of  the  parameters  Tv,  8V  as 

2  p  D  L 

given  by  ^  T  cr  or  ^  8V  \s  d-  'j 

Thus  Fig*  22  indicates  that  the  Bagnold's  theory  can  be 
applied  to  the  behaviour  of  practical  suspensions  in 
capillary  viscometer  with  suitable  modifications*  Also, 
it  indicates  that  the  dimensionless  terms  obtained  by 
dimensional  analysis  adequately  represent  the  phenomenon* 

It  is  realised  that  in  the  present  Investigation 
the  size  of  the  sand  d,  the  density  of  sand  ^  s,  and  the 
viscosity  of  water  ,  were  constant,  and  hence  the  product 
?,a2  was  a  constant*  However  the  encouraging  correlation 
of  Fig*  22  and  its  similarity  with  Bagnold’s  experimentally 
derived  curve  (Fig*  6)  indicate  that  the  dimensionless  terms 
used  are  meaningful* 


Figs.  17,  18,  indicate  that  while  the  behaviour 
of  the  suspension  in  the  concentration  range  0.13  to  0*30 
follow  a  particular  trend  there  is  a  sudden  change  in  the 
characteristics  between  C  “  0*30  and  0*35*  Even  though  the 
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region  C  -  0*13  to  0*30  was  the  region  of  interest  in 
the  previous  analysis  the  experimental  points  with 
C>0*30  are  all  Included  in  the  analysis  and  the  plots* 

It  i 8  clear  from  figs*  17  through  22  that  there  exists  a 
separate  regime  of  behaviour  for  C>0*30  but  the  exper- 
imental  data  are  insufficient  to  predict  the  behaviour 
In  this  region* 

Also,  the  existence  of  a  critical  concentration 
between  C  ■  0*30  and  0*35  is  clear  from  the  above  plots* 
This  phenomenon  of  critical  concentration  has  been  observed 
by  several  investigators*  The  various  observations  and 
the  explanations  advanced  for  this  behaviour  have  been 
reviewed  in  sec*  2*4*  The  observed  value  of  Ccr  between 
0*30  and  0*35  is  within  the  range  of  values  of  0*30  and 
o*50  at  which  this  phenomenon  has  been  observed  by 
various  investigators* 


' 

4  1  * 

' 


CHAPTER  V 


CONCLUSIONS  AND  RECOMMENDATIONS 
CONCLUSIONS 

The  capillary  viscometer  developed  in  the 
laboratory  is  well  suited  for  the  analysis  of  the 
rheological  characteristics  of  suspensions*  The  tests 
on  the  sand-water  suspension  show  that  systematic  test¬ 
ing  can  be  carried  out  over  a  very  wide  range  of  shear 
rates  in  such  equipment  to  get  useful  and  sufficiently 
accurate  data* 

The  review  of  the  literature  on  flow  characteristics 
of  suspensions  indicate  that  the  available  information  is 
insufficient  to  predict  the  performance  of  a  suspension  in 
a  practical  case*  Also  in  most  of  the  theoretical  and  empirical 
works  there  is  a  lack  of  full  appreciation  of  some  of  the 
factors  influencing  the  behaviour  of  the  suspension*  The 
Bagnold's  concept  of  grain  shear  stress  is  a  new  approach 
and  seems  to  have  considerable  potential  of  application* 

The  tests  on  sand-water  suspension  indicate  that 
In  the  range  of  shear  rates  2300  -  6300  sec~*  the  suspension 
is  non-Newtonian  dilatant  at  all  concentrations  greater  than 
0.13*Using  Metzner#s  method  of  describing  the  capillary 
viscometer  data.  It  is  found  that  the  suspension  behaves  as 
a  powerlaw  fluid  with  n  »  4/3  for  values  of  C  between  0*13 
to  0*30* 

The  results  of  the  dimensional  analysis  when 
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applied  to  the  test  data  indicate  that  the  dimensionless 

2  T  ?  d'2' 

number  Gz,  -  — -  has  a  unique  functional  relation- 

'T  A 

ship  with  the  term  (  %v  d  >  for  a  concentration  range 

X  <  3*5*  By  considering  the  similarity  of  this  functional 

2 

relationship  and  that  between  and  M  as  found  by 

Bagnold,  a  critical  number  G)|  cr  can  be  envisaged  to  mark 
the  transition  from  Newtonian  to  non-Newtonian  flow 
regime  of  the  suspension* 

A  critical  concentration  at  which  an  abrupt 
change  in  the  characteristics  of  the  suspension  take 
place  is  noticed  between  C  *  0*30  and  0*35* 

It  is  recommended  that: 

(a)  a  systematic  study  be  carried  out  by  varying  the 
parameters  over  a  wider  range  to  prove  conclusively 
the  validity  a4kd  limitations  ot  the  above  conclusions* 

(b)  a  similar  study  be  made  for  turbulent  flow  through 


vertical  tubes* 


. 
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APPENDIX  A 


A-l 


TABLE:  A-l 

CALIBRATION  OF  CAPILLARY  TUBE:  (A  -  M) 


TUBE:  A, 

LENGTH  = 

122  cm.s 

DIA:  2mm 

OIL:  M, 

TEMP.  « 

72°F 

VISCOSITY : 

U.  I4I6  c  .p  . 

SP.GR. :0.85U7 

EXPT.  NO. 

1 

1 

j 

VOL. 

c.c. 

TIME 

secs 

PRESS 

DIFF. 

Ctr\. 

A'p 

u  -  1 

T 

UJtri 

•Dyines/  Cm1 

8V 

D 

-1 

sec 

1 

1*3.0 

65.3 

11*1.0 

U8.5 

8I4.O 

2 

1*7.9 

71.6 

11*1.0 

1*8.5 

870 

3 

1*2.0 

29.1 

295.0 

101.5 

1835 

j  h 

39.0 

27.2 

295.0 

101.5 

1830 

5 

U6.5 

28.3 

3U0.0 

115.0 

2100 

6 

1*6.0 

29.2 

317.0 

109.0 

2000 

7 

1*5.0 

25.1 

365.0 

125. 5 

2280 

8 

1*2.0 

20.8 

1*01.0 

138.0 

2580 

9 

U1.0 

17.9 

1*59.0 

158.0 

2930 

10 

ia. 5 

16.0 

508.5 

175.0 

3310 

11 

1*2.0 

12.8 

670.0 

230.0 

I4I8O 

12 

1*6.0 

11.3 

815.0 

280.5 

5200 

A- 2 


TABLE:  A  -  2 

CALIBRATION  OF  CAPILLARY  TUBE:  (A  -  N) 

DIA  «  2  mm. 

VISCOSITY  * 
9.017  c.p. 

SP.  GR.  * 
0.851 


EXFT. 

NO, 

VOL. 

c.c. 

TIME 

secs 

PRESS 

DIFF. 

Hm  | 

T 

1  Cvrp- 

8V 

D 

-1 

sec 

1 

h2.0 

123.li 

136.5 

I16. 7 

U33 

2 

39.0 

52. 1 

30U.0 

lOli.O 

955 

3 

U7.5 

U6.5 

1*22.0 

1I1U.0 

1300 

Ui.o 

36.0 

503.0 

172.0 

1555 

5 

lo..  8 

28.3 

585.0 

200.0 

1885 

6 

U2.5 

2U.6 

70U.0 

210.0 

2200 

7 

10.8 

20.  k 

825.0 

282.0 

2620 

8 

1*6.8 

19.6 

957.0 

328.0 

3010 

TUBE:  A  LENGTH  -  122  cm. 

OIL:  N  TEMP.  «  71°F 


• 

A- 3 


TABLE?  A  -  3 

CALIBRATION  OF  CAPILLARY  TUBE:  (B  -  M) 

DIA:  2.7  mm, 

VISCOSITY: 
1+.1+16  c.p. 

SP.  GR.  0.851+7 


EXPT. 

NO. 

VOL, 

c.c. 

TIME 

secs 

PRESS 

DIFF. 

H  cm 

m 

T 

v  WTr* 

"Dynes  /  cm1' 

8V 

D 

-1 

sec 

1 

U6.5 

17.8 

139.0 

66.8 

1350 

2 

1+6.5 

11+.1+ 

170.5 

81.8 

1670 

3 

1+5.5 

11.0 

218.0 

iol*.5 

211*0 

1+ 

1*3.0 

8.9 

261.0 

125.0 

2500 

5 

1+6.5 

8.1+ 

290.0 

139.0 

2860 

6 

l+l+o  0 

5.6 

1*21.0 

202.0 

1+060 

7 

1*9.5 

18.8 

138.5 

65.5 

1360 

8 

1*5.0 

17.0 

139.1 

67.0 

1370 

9 

1+6.5 

11.2 

219.5 

105.0 

2150 

10 

l*5.o 

8.1 

295.0 

11+2.5 

2870 

11 

1*1*.  5 

6.1+ 

366.0 

175.5 

3590 

12 

85.0 

7.8 

1+58.0 

220.0 

1*530 

13 

67.0 

7.1+ 

1*73.0 

227.0 

1+690 

TUBE:  B  LENGTH:  118  cm. 

OIL:  M  TEMP.  «  72°F 


' 


A-4 


TABLE.*  A  -  k 

CALIBRATION  OF  CAPILLAR!  TUBE?  (B  -  N) 

DIA?  2.7  mm. 

VISCOSIT I ? 
80872  c.p. 

SP.  GR.  0.8508 


EXPT, 

NO. 

VOL. 

C*  Co 

TIMS 

secs 

PRESS 

DIFF. 

Hjjj  0, 

T 

V  COYt*| 

l^yne^  1  c*Y*a* 

8V 

D 

-1 

sec 

1 

U5.5 

31*.  7 

11*1.5 

68.0 

679 

2 

U3.0 

17.2 

261.0 

125.2 

1300 

3 

1*6.5 

xii.  3 

352.0 

169.0 

1680 

h 

l*o.o 

9.9 

1*33.0 

208.0 

2090 

5 

38.5 

8.5 

1*70.0 

226.0 

2350 

6 

1*3.5 

8.2 

559.0 

268.5 

271*0 

7 

1*6.5 

9.5 

£20,0 

250,0 

251*0 

8 

1*9.0 

11*.  1* 

358.0 

172,0 

1760 

9 

1*8.0 

22,6 

222.5 

107.0 

1100 

10 

1*8. o 

35.3 

11*0.5 

67.5 

705 

_ 1 

TUBE?  B  LENGTH?  118  cm. 

OIL?  N  TEMP.  -  72°F 
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TABLE:  A-5 

SAND-WATER  SUSPENSION  TEST  DATA  -  A. 


SERIES:  A  D  «  0,20  mm*  oc  ®  0«82  L  ■  122  cm, 

-  9©1  x  ^  ,  Oi 

8V  «  -S  Y-  i 

IT  IiH  "  ^ 


EXPT, 

NO, 

VOL. 

V  c.  c. 

TIME 
t  sec. 

DIFF. 

PRESS 

HEAD 

VArrv 

CONC. 

C. 

SP,  GR. 
SG 

T 

X>'-/»V2.S  j  Cm1 

8V 

TT 

sec"^ 

A-l 

1*3-5 

18.6 

166 

0.125 

1.205 

65.9 

2978 

A- 2 

1*2.2 

11*.  6 

218 

0.129 

1.212 

87.1 

3680 

A- 3 

kk.o 

12.0 

310 

0.129 

1.212 

123.9 

1*669 

A-U 

Uo.5 

12.3 

262 

0,13k 

1.220 

105.1* 

1*192 

A-5 

UU.5 

10.55 

381 

0.131* 

1.220 

153.2 

5371 

A- 6 

1*5-0 

22.0 

11*2 

0.159 

1,260 

59.0 

2600 

A- 7 

UU.5 

18.3 

185 

0,159 

1.260 

76.8 

3096 

A- 8 

1*3-9 

13.63 

256 

0.171* 

1.275 

107.6 

1*101 

A- 9 

35-0 

11.80 

228 

0.201 

1.330 

100,0 

3777 

A- 10 

1*3-0 

17.20 

216 

0.216 

1-355 

96.2 

3183 

A- 11 

37-5 

19.70 

16U 

0.226 

1-370 

741 

21*21* 

A- 12 

1*3.5 

20.10 

175 

0.250 

l.l*lo 

81.3 

2756 

A-13 

1*3.0 

9.90 

1*90.5 

0.271 

1. 1*1*5 

233.7 

5530 

A-1U 

1*1.5 

10.10 

U88 

0.277 

1-1*55 

23l*.l 

5232 

a-15 

i 

1*1.6 

8.65 

559 

0.278 

1.1*55 

268.1 

6123 

A-6 


TABLE;  A~5»  continued 


EXPT, 

NO. 

VOL, 

V  o.  c. 

TIME 
t  sec. 

DIFF. 

PRESS 

HEAD 

H*i  Cm 

CONG. 

C 

SP.  OR. 
SG 

T 

v  00 

8V 

IT 

sec“^ 

A- 16 

UUoO 

8.80 

67U 

0.287 

1.U70 

326.6 

6366 

A-17 

1*7.5 

1*8,60 

178 

0.299 

1.1*90 

87.1 

121*1* 

A- 18 

1*1.  U 

8.90 

590 

0.302 

1.1(95 

290.8 

5923 

A- 19 

111*.  8 

10.80 

52U 

0.305 

1.500 

259.1 

5282 

A- 20 

U2.2 

ll*.80 

310 

0.310 

1.510 

1511.3 

3630 

A- 21 

Uo.o 

19.8 

291 

0.317 

1.520 

11*5.8 

2572 

A- 22 

1(1.6 

8.3 

71(1 

0.323 

1.530 

373.8 

6382 

|  A-23 

35.0 

103.1* 

13U 

0.326 

1.535 

67.8 

1(71 

k-2h 

31.0 

58.6 

208 

0.332 

1.5U5 

105.9 

67U 

A-25 

32.0 

102.8 

13U 

0.335 

1.550 

68.5 

396 

A-26 

33.0 

103.2 

13U 

0.335 

1.550 

68.5 

U07 

A-27 

3U.0 

38 .3 

382 

0.335 

1.550 

195.2 

1130 

A-28 

1(3.0 

51.7 

336 

0.338 

1.555 

172.3 

1059 

A-29 

1*1.0 

58.6 

3U8 

0.3U1 

1.560 

179.0 

892 

A“30 

L- 

1(3.5 

1(7.2 

399 

O.3I4I 

1.560 

205.2 

1173 

A- 7 


TABLE  A  -  continued 


EXPT. 

NO. 

VOL. 

V  c.c. 

TIME 
t  sec. 

DIFF. 

PRESS 

HEAD 

Hw\  Cvri 

CONG. 

C 

SP.  GR. 
SG 

Xo 

'D'yv'ie.t,  j  cm1 

8V 

IT 

S9C"1 

A- 31 

29.1 

100.9 

13U 

0.31*2 

1.562 

69.0 

367 

A- 32 

1*3.0 

37.8 

1*65 

0.31*2 

1.562 

239.2 

114*8 

A- 33 

39.2 

65.3 

239 

0. 3J4.U 

1.565 

123.3 

76U 

A -3k 

1*1.8 

51.8 

3U7 

0.31*5 

1.565 

179.0 

1027 

A-  35 

38.0 

U8.7 

316 

0.31*9 

1.572 

163.8 

991* 

A- 36 

35.5 

92.7 

150 

0.353 

1.580 

78.1 

488 

A- 37 

311.2 

71.5 

186 

0.353 

1.580 

96.8 

609 

A- 38 

32. 5 

95.6 

13U 

0.357 

1.585 

70.0 

1*33 

1 

TABLE:  A  -  6 


SAND- "WATER  SUSPENSION  TEST  DATA  -  B 


SERIES:  B  D  °  0,27  cm*  o4  «  0.90  L  «  11.8  cm. 


EXPT. 

NO. 

VOL. 

Vc.c. 

time 

t  sec. 

DIFF. 

PRESS 

HEAD 

Hr**  Cw\ 

l 

GONGN. 

G 

SP.  GR. 
SG 

T 

V 

8V 

D 

B-l 

5o.o 

U.  30 

310 

.136 

1.220 

190.5 

6017 

B-2 

105.0 

8.10 

371 

.151 

1.250 

23a.  2 

6708 

B-3 

67.8 

6.80 

21*1 

.161* 

1.270 

15a.  0 

5160 

B-U 

71.2 

6.20 

277 

.161* 

1.270 

176.0 

59a3 

B-5 

5i.o 

7.20 

160 

.166 

1.275 

102.6 

3666 

B-6 

69*7 

5.90 

338 

.182 

1.290 

220.2 

6111* 

B-7 

65.7 

6. 14.0 

259 

.181* 

1.300 

170.0 

5312 

B-8 

69.7 

7.30 

251 

.187 

1.305 

165.  a 

a9ai 

B-9 

53.6 

7.30 

170 

.189 

1.310 

112.2 

3800 

B-10 

81.3 

9.30 

203 

.192 

1.325 

135.8 

h$2h 

B-ll 

60.0 

6.70 

228 

.191* 

1.315 

151.1* 

1*631* 

B-12 

58.3 

7.10 

206 

.196 

1.320 

137.3 

a2a9 

B-13 

80.7 

7.1*0 

306 

.196 

1.320 

20i*.0 

561*1* 

B-ll* 

70.5 

5.25 

1*15 

.199 

1.327 

278.0 

6955 

B-15 

1*6.2 

7.2 

137 

.204 

1.335 

92.7 

3321 

1 

j! 

1 

-J 

A- 9 


TABLE;  A  -  6,  continued 


EXPT , 

NO. 

VOL. 

V  c.c. 

TIME 
t  sec. 

DIFF. 

PRESS 

HEAD 

Hrrx  Cvn 

GONCN. 

C 

SP.  GR. 
SG 

t 

Vuo 

8v 

t r 

B-16 

78.7 

13.65 

139.5 

.201* 

1.335 

93.7 

2981* 

B-17 

UU*o 

7.60 

mo 

.201* 

1.335 

9l*.l* 

2996 

B-18 

81.3 

12.0 

167 

.201* 

1.335 

112.2 

3506 

B-19 

1*1*.  1* 

3.50 

380 

.222 

1.365 

262.0 

6565 

B-20 

61.3 

11.00 

137 

.235 

1.380 

95.5 

2881* 

B-21 

69.5 

10.80 

158.5 

.238 

1.380 

110.5 

3330 

B-22 

1*0.5 

5.oo 

231 

.250 

1.390 

162.2 

1*192 

B-23 

h6.5 

6.30 

179 

.257 

I.I4IO 

127.5 

3820 

B-2i* 

61.0 

8.70 

193 

.257 

1.1*20 

128.7 

3628 

B-25 

82.7 

9.20 

2U1 

.257 

1.1*20 

172.8 

1*650 

B-26 

81.5 

7.95 

295 

.257 

1.1*20 

211.6 

530.5 

B-27 

1*7.0 

8.50 

11*2 

.259 

1.1*25 

102.2 

2862 

B-28 

1*3.0 

6.00 

198 

.262 

1.1*30 

1 1*3*0 

3708 

B-29 

l*o.o 

6.50 

163 

.263 

1.1*31 

115*6 

3185 

B-30 

1*3.0 

5.8o 

211* 

.268 

1. 1*1*0 

155.6 

3837 

B-31 

1*3.0 

l*.5o 

312 

,280 

1.1*60 

230.0 

1*91*5 

B-32 

l*5.o 

5.30 

256 

.283 

1.1*65 

189.1* 

1*391* 

B-33 

58.0 

5.io 

388 

.283 

1.1*65 

286.1 

5885 

B-3iU 

1*3.0 

7.1*0 

170 

.292 

l.l*8o 

127.1 

3007 

• 
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APPEflMX  -  £ 

RELATIONSHIP  BETWEEN  X  AND  C, 

Consider  an  array  of  spheres  of  diameter  d  in  a 
known  type  of  piling,  packed  such  that  minimum  amount  of 
voids  corresponding  to  the  chosen  piling  are  produced* 

Let  V=>  ■  Volume  of  Solids 

V  ■  Total  Volume 

V5 

C  -  —X  -  Volume  Concentration  of 

*  V  Solids. 

In  this  packing,  since  the  amount  of  voids  is  minimum,  the 

free  distance  between  adjacent  spheres  is  zero  and  hence  by 

definition  (Sec*  2*3)  X-  °°  *  Hence  C*  is  the  concentration 

I 

at  X  „  00  * 

Now  let  the  distance  between  the  adjacent  spheres 
be  increased  to  bd,  keeping  the  piling  arrangement  the  same* 

In  a  total  volume  of  V , the  amount  of  solids  present 

V 

is  now  — 1  ,  The  resulting  concentration  of  solids  is 

D 

C  -  -  '  -%L  - (B-l) 

b3V  b5 

But  from  Fig#  5,  bd  -  a  +  d 

b  «■  A  +  1 
d 

-  1/x  +1  - (B-2) 

c  -  —  „ 

c-k  +  |)’ 

or  A  -  - - - - 

(f)1/3  -1  . (B-3) 


■-  -  •  .  ,  5  -  t  i 


In  the  above  Equation  B-3  C*  is  the  maximum 
possible  concentration  of  spheres  in  a  given  type  of  piling. 
The  value  of  C*  fox  various  types  of  pilings  has  been 
calculated  by  Slichter  (25).  It  is  found  that  the  least 
valve  of  C*  is  0*524  for  cubic  piling  and  the  greatest  value 
of  C*  ■  0*74  for  rhombohedral  piling*  Bagnold  (20)  has 
found  C*  ■»  0*65  for  natural  uniformly  rounded  sand  grains 
under  dry  compaction* 

Fig*  23  shows  the  relationship  between  X  and  C 
Calculated  from  equation  B-3  with  C*  -  0*65*  This  Fig.  is 
used  in  the  calculation  of  X  corresponding  to  various 
values  of  C  in  Tables  A-7  and  A-8* 
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